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LXXXIX. On the Theory of Beta-Radioactivity 


By A. M. Suira 
Department of Natural Philosophy, University of Aberdeen* 


[Received April 15, 1952] 


ABSTRACT 


A comparison of some experimental beta spectra with theoretical 
predictions when the interaction is taken as a linear combination of the 
five usual forms has been made, and on this basis the conclusion is that 
the interaction can only be a mixture of the vector and tensor forms. 
. For two disintegrations involving light nuclei the arbitrariness in the 
theoretical predictions has been reduced by assuming a shell structure 
of nuclei and evaluating the nuclear matrix elements involved. From 
this analysis it is concluded that the interaction responsible for beta 
decay is a mixture of approximately two parts of vector form to one of 
tensor. 


$1. INTRODUCTION 


In the Fermi theory of beta radioactivity the interaction between light 
and heavy particle fields is expressed in terms of Dirac’s four-component 
relativistic wave functions for the particles involved in the decay process, 
in such a way that the interaction is invariant under the Lorentz group of 
transformations. There are five such expressions, called the scalar, 
vector, tensor, axial vector and pseudoscalar forms of interaction (S, V, 
T, A and P in the usual notation). The theory of allowed and forbidden 
transitions has been developed, for each of these five different forms, by 
Konopinski and Uhlenbeck (1941). For allowed transitions all forms of 
interaction predict the same distribution of energy of the emitted beta 
radiation ; for forbidden transitions the results are expressed in terms of 
correction factors by which the allowed distribution must be multiplied. 
Each correction factor contains a few nuclear matrix elements, each 
multiplied by a coefficient which is dependent on the energy of the emitted 
beta particle, and comparison of theory with experiment is rendered 
_ difficult by lack of knowledge concerning the matrix elements. However, 
the selection rules which the matrix elements impose on the spin and 
parity changes in the disintegration do limit the correction factors in a 
few cases to only one matrix element along with its energy dependent 
coefficient. For example, in the first forbidden correction factor for the 
tensor or axial vector interactions, the matrix element B,; allows a spin 
change of 2 units and a parity change, whilst the other matrix elements in 
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this correction factor are zero under these circumstances. Now, several 
spectra have been measured experimentally such that when the correction 
factor (B,; shape) (the energy dependent coefficient of 2'|B,,; |?) is applied 


i,j 
to the conventional Fermi plot, a straight line is obtained. For example 
there are the spectra of ®S1, °°Sr, °Y and *Y (Langer and Price 1949). 
Moreover, the Fermi plot of the spectrum of !°Be gives a straight line 
when corrected for by (Sj, shape)—the second forbidden correction 
factor for tensor or axial vector interactions when a spin change of 3 
units and no parity change are involved (Fulbright and Milton 1949) ; 
and so also with #°K when corrected for by the so-called C factor, the 
coefficient of the matrix element allowing a spin change of 4 units and a 
parity change in the third forbidden tensor or axial vector correction 
factors, thereby giving a unique energy dependence (Alburger 1950). 
Thus the interaction must be, partly at least, of the tensor or axial vector 
form. 

§2. A LINEAR COMBINATION OF INVARIANTS 


There is no theoretical reason why one form of interaction should be 
preferred to any other, and so the additions to the correction factors, 
which arise when a linear combination of the five invariants is considered, 
have been worked out by the author and are quoted elsewhere (Smith 
1951). A linear combination of interactions has been suggested by 
Fierz (1936) who investigated the effect such a combination would have on 
the allowed distribution. He showed that if the scalar interaction is 
combined with the vector, or the tensor with the axial vector, the allowed 
distribution must be multiplied by (1-++c/w) where c is a constant. Since 
many beta spectra conform to the allowed distribution as predicted by 
Fermi, such combinations must be impossible. 

We take then as the invariant 


with the understanding that 
AgAy= ApA4g= 0. 


The constants A are taken to be real. 

It is worth noting that those explanations of spectral shapes which 
have already been established would not be affected by such a combination 
of interactions. (There has been a recent exception to this—it is possible 
to explain the shape of the spectrum of *°Cl using a mixture of coefficients 
of those matrix elements in the second forbidden tensor correction factor 
which allow the known spin change of 2 units. This will be dealt with 
more fully in §5.) Allowed spectra are not affected. Spectra which 
are explained by (B;; shape) would not be affected, since there is no matrix 
element in the first forbidden correction factors of either the scalar or 
vector interactions which allows a spin change of 2 units. The same is 
true of spectra explained by (S,;, shape) and by the C factor. It is therefore 
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reasonable to try a combination of interactions with a view to explaining 


the spectra of nuclei which do not conform to the distributions for single 
interactions. 


§3. COMPARISON WITH EXPERIMENT I 


Comparison of theory with experiment is made difficult by the lack of 
reliable spectra whose shapes have not hitherto been explained. It has 
been realized in recent years that a considerable amount of distortion can 
occur at the low energy end of a beta spectrum unless the emitting layer of 
radioactive material is extremely thin (Albert and Wu 1948). However, 
the literature has been searched for unexplained spectra which seem to 
be reliable, and six disintegrations have been selected in which (a) the 
source used was reasonably thin, (6) the ft value suggests a second 
forbidden transition. Since a spin change of 3 units would give a spectrum 
conforming to (S,;;, shape), and a spin change of 1 unit would give an 
allowed transition, we assume that a spin change of 2 units and no change 
in parity are involved. These assumptions are justified in some cases by 
predictions of the shell model of nuclei, but should it be impossible to 
explain any spectrum under these assumed conditions, it need not be 
taken as a breakdown of the theory. Furthermore, this procedure will 
not indicate whether the interaction contains the pseudoscalar form, 
since it gives no contribution to second forbidden transitions with a spin 
change of two units. 

The shape which the correction factor should have is determined by 
finding the factor by which the conventional Fermi plot must be corrected 
in order to give a straight line. In this section we examine whether it is 
possible to combine the coefficients of the matrix elements occurring in the 
complete correction factors so as to give the experimentally determined 
correction factors, without making any assumptions concerning the 
magnitudes of the matrix elements or the proportions in which the 
interactions are mixed. 

The fitting of a combination of so many theoretical curves on to one 
experimentally determined curve is done by the method of least squares. 
We have to fit a combination of f,, f., fz, etc., on to an experimentally 
determined curve f and this is done by choosing aj, a3, a3 etc., so that 


Aah, + efe,+Gfs, Sire aly 


is a minimum, the summation being over points on the curves. This 
gives rise to a system of linear equations which can be solved for the 
constants @,, @y, a3, etc. However, some of the functions /;, f,, fs etc. 
usually have similar shapes, and consequently there will be approximate 
linear relationships between them. This will make the matrix of the set of 
linear equations nearly singular, and the number of functions /;, f2, fs - - - 
must be reduced so that only independent ones are used in the fitting ; 
otherwise it might be possible to get a good fit which would not be allow- 
able, since it might be using the small errors in the linear relationships as 


3R2 
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independent functions. Moreover, because of this singularity, we will not 
be able to deduce anything about ratios of matrix elements or the pro- 
portions in which interactions are mixed; we may be able, however, to 
eliminate certain combinations of interactions as being quite incapable of 
explaining the spectra considered. 

Evaluation of the matrix element coefficients in the correction factors, 
and the subsequent curve-fitting processes require a very large amount of 
computational work. For this purpose we require to express the ratios 
F,/F and F,/F* which appear in the correction factors in forms suitable 
for computation. 

These are 

Fy, (4!)? F? (28-+-1) 


F > 48,41 (2pp)"—8—? [(S.—1*° +9") 


02Z2 a2 Z2y? & l 
eS, oh : 
ss E 2 oa 2 n=1 n(n*+-y") 
Cet BCE es 5 ee HC aera ee 
Fike WE (28,1) (2pp) [(Sg—1)?+-y¥*]L(S2—2)?+y"] 
| 2427 2a2Z2y? a) 1 
x} i- 


BOT ary 
Here the symbols have their usual meaning, C is Gauss’s number=0-577 
and y=aZW/p. The nuclear radius p is given, in relativistic units, by 
p=3-67 x 10-3x Al8 

where A is the mass number. 

For each disintegration the best fits for the four possible combinations, 
S or V with T or A have been worked out. In addition, the best fit for a 
pure tensor interaction has been found, this being the only interaction 
which by itself gives promise of explaining all spectra. 


Phosphorus 

The PS disintegration has been reported as giving a straight line 
Fermi plot by Siegbahn (1946) and by Langer and Price (1949), but a 
more recent determination by Agnew (1950) gives a Fermi plot which is 
slightly concave towards the energy axis and has a sudden rise at low 
energies. A later measurement by Warshaw, Chen, and Appleton (1950) 
also gives this sudden rise at low energies but gives a straight line Fermi 
plot above this. The ft value, based on Agnew’s determination of the 
end-point energy, is 8-6 x 107, which classifies the transition as second 
forbidden, and the conditions 4J= 2, no change of parity, are quite reason- 
able on the shell model of nuclei. Since Agnew used a very thin source 
(<0-04 mg/cm*) and since the departure from a straight line Fermi plot, 
apart from the sudden rise in low energies, is in the opposite direction to 
what one would expect from source thickness distortion, his results have 
been selected for comparison with second forbidden correction factors. On 


ORS Se SANs Madara Pi ce SL ieee 


* Konopinski and Uhlenbeck, loc. cit., Smith, loc. cit. 
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the assumption that the sudden rise at low energies does not imply 
complexity of the spectrum, it has been found impossible to fit the experi- 
mental data with any of the second forbidden correction factors, and we 
must conclude that °2P does not fall in this class of transitions. It is 
probable that this transition is really allowed, in accordance with all the 
other experimenters, with a spin change of 1 unit. The ft value would then 
be high because of some inner selection rule which is making the allowed 
matrix element fo small (See Pursey 1951 a).* 


Technetium 
The *Tc+*Ru spectrum has been measured recently by Taimuty 
(1951) using a fairly thin source (0-17—1-3 mg/em?). He found that the 
Fermi plot was reasonably good when corrected for by (B;; shape), which 
would make the transition first forbidden with a spin change of 2 and a 


Fig. 1 


1-0 - 1:2 Ww }-4 16 
°°Tc. Best fit for combination S7’. 


change of parity. The ft value, however, is 2-310! which indicates a 
higher order of forbiddenness, and the shell model suggests a spin change 
of 2 units with no parity change. Attempts to fit the spectrum with 
second forbidden correction factors show that it is possible to do so with 
the tensor-interaction alone or with any of the four possible combinations 
of interactions. A typical fit is shown in fig. 1 where the two dotted 
curves represent the experimentally determined correction factor for the 
experimental limits put on the end-point energy, and the solid curve is the 
best fit obtained from the theoretical curves by the method of least 
squares. 
‘Tin 

The disintegration ”5Sn+>"5Sb has been investigated by Hayward 

(1950) using a source of thickness estimated at 0-4 mg/cm. The ft value 


Oo Ee ey ci ee ees eee eee 2 

* The work of Jensen et al. (1952) has just come to the author’s notice. They 
find the spectrum of *?P to be allowed, and suggest that the sudden rise at low 
energies is due to the presence of the **P isotope. 
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is 3:7 108, which puts the transition in the second forbidden class. The 
decay is obviously complex. Hayward has corrected the usual Fermi 
plot with (B,; shape) and claims the resulting plot to be a straight line ; 
however, examination of it reveals that it has a definite convexity towards 
the energy axis. If we assume a spin change of 2 units and no parity 
change, which is quite reasonable on the shell model, we find that it is 
possible to explain this spectrum with the tensor interaction alone or with 
any of the four combinations. 
Thulium 

17”Tm decays by beta emission to 1Yb and this disintegration has also 
been studied by Agnew who reports it along with *P, and the Fermi plot 
he obtains is slightly concave towards the energy axis. The spectrum has 
been measured by a few other experimenters ; Saxon and Richards (1949 a) 
found that it was of the allowed form ; Fraser (1949) considered the decay 
to be complex, being made up of two transitions of the allowed form ; 
Graham and Tomlin (1949) found it to be allowed. The ft value is 7 x 108, 
and since Agnew used a very thin source (<0-03 mg/cm*), his results 
have been compared with second forbidden correction factors. None of 
the interactions gives a satisfactory fit. We must conclude that the spin 
and parity assumptions are not valid for this disintegration, and that it 
must either be allowed or first forbidden, in accordance with the measure- 
ments of Saxon and Richards and of Graham and Tomlin. The deviation 
from the allowed form in Agnew’s measurement is very similar to that in 
his determination of the spectrum of **P. 


Thallium 
The beta spectrum of 2°T]+2%4Pb has been measured by Saxon and 
Richards (1949 b) using a source of thickness 0-1 mg/em?. They have 
corrected the Fermi plot obtained by the (B,; shape), and the resulting 
plot gives a straight line at the high energy end but becomes convex 
towards the energy axis at lower energies. This deviation, they consider, 
cannot be attributed to residual scattering. The ft value is 5x 10°, 
putting the transition in the second forbidden class. Comparison of the 
theoretical correction factors with the experimental form shows that it 
is possible to explain this spectrum well by tensor interaction alone, by 

ST’, VT or VA, but not by SA (see fig. 2). 


Radium E 

Konopinski and Uhlenbeck (1941) have explained the spectrum of 
RaE using the second forbidden correction factor for the tensor interaction. 
The experimental results were by Flammersfeld and also by Langer and 
Whitaker, but the fit obtained does not agree so well with the experimental 
work of Neary. Neary’s results have since been confirmed by Langer 
and Price (1949) and the measurements of Langer and Price have been 
used here for comparison with theory. The ft value is 9-8 x 107, which is 
consistent with a second forbidden transition. It has been found possible 


/ 
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to explain this spectrum with interactions SA and V7 only, the others 
giving a definite divergence from the experimental curve. Figure 3 (a) 
shows the best fit for the tensor interaction and for ST ; the best fit for VA 
is sumilar to that for the tensor interaction. The best fits for V7 and SA 
are almost identical and are shown in fig. 3 (b). 


Fig. 2 


1-0 Sumber Grpea:0 25 


20471]. Curve (1) is best fit for SA; curve (2) is best fit for ST. Best fits for 7, 
VT, VA are similar to (2). 


6) 
1:0 5 W 2:0 2:5 3-0 ie) is Ww 2-0 2:5 0) 


(a) (b) 
Rak. Best fit for 7 or ST. Rak. Best fitfor V7 or SA. 


The results of this analysis are summarized in table 1, the Yes or No 
indicating whether a spectrum can or cannot be explained by a particular 
interaction. Excluding **P and 'T'm from consideration as second 
forbidden transitions, we conclude that if our assumptions concerning the 
spin and parity changes involved in the other disintegrations are correct, 
then the interaction can only be a mixture of the vector and tensor forms. 

One other point emerges from these considerations—the ft value should 
be looked upon as setting an upper limit to the order of forbiddenness of a 
transition rather than giving the actual order. The ft values of *P and 
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170Tm suggest that they are second forbidden transitions whereas they 
appear to be allowed or at most first forbidden. In the same way the ft 
value for “Ci (~7 x 108) classifies its decay as second forbidden although 
its spectrum has recently been found to be of the allowed form (Feldman. 
and Wu 1949). 


Table 1 


§4. THe Nuctear Matrix ELEMENTS 
The nuclear matrix elements are of the form 


z i V*0,Q,Udr 
k 


where U and V are the wave functions of the initial and final nuclei 
respectively. @), is an operator operating on the isotopic spin of the kth 
nucleon in the initial nucleus in such a way that it changes the character 
state of a nucleon from a neutron state to a proton state. O is an operator 
having the forms given by Konopinski and Uhlenbeck (1941). In their 
paper they have dropped the operator f, but its retention is essential since 
it changes the sign of matrix elements like R,;, 7';; which combine large 
components of the initial wave function with large components of the final 
one, and the sign is of great importance in cross correction factors ; 
moreover, since 8 does not commute with «, its presence can change 
completely the values of matrix elements such as A,,. We denote by 
Af,, Tf, ete., those matrix elements which are the same as A go Ty ete., but 
have the operator 8 included in the integrand. 

In recent years, a considerable amount of evidence has been advanced in 
favour of the shell structure of nuclei, and this model is used here to work 
out the matrix elements more exactly. The nucleons are considered as 
individual particles moving in an average field due to all the others, and so 
the nuclear wave function will be a sum of products of the nucleon 
individual particle wave functions, so arranged as to be an eigenfunction of 
the angular momentum J and of the magnetic quantum number M. We 
are interested in the matrix elements for a transition from an initial nucleus 
with angular momentum J; and magnetic quantum number M ; to a final 
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nucleus with angular momentum J; and any magnetic quantum number 
M;, and so we must sum the transition probabilities over all permissible 
values of M,. Thus, if P(J;, M;; J;, M rs) be the transition probability 
from an initial state J;, M; to a final state J;, M;, and Q be the total 
transition probability to all values of M;, then 


oF 
QI, Mi35,)= 2 PU, M;; Jy, My) ae = eae: 
My=—Jf 
there being (2J,+1) terms in this summation. Now these transition 
probabilities have been expressed in terms of the matrix elements of 
irreducible tensors, along with their energy dependent coefficients, and 
Wigner (1931) has shown that if bee pe ae be the matrix element of 
the p-component of an irreducible tensor for the transition from an 
initial state J;, M; to a final state J,, M;, then 
1 
@ pee o 2 ¢ 
Patera - Spare Pag Mis Tp My Peace) 
The right-hand side of (2) is independent of M,, and so Q in (1) can be 
written as Q(J;; J;). In other words, the transition probability is inde- 
pendent of our choice of axes, as we would expect. Furthermore, since 
the operators in the matrix elements in beta theory are Hermitian, we 
obtain from (2) 
(2J,+1) 2 [T5,, my; 3,00, P= (27 +1) DH [T5,,0;; sym, P- 
e, M; Q, M; 
From the derivation of Wigner’s result, it is obvious that this applies also 
to scalar products of matrix elements such as 2 A,,* 7’. We have, 
therefore, i 
Pe Oy TN). Ig. J i, «, (8) 
Thus in working out matrix elements we are at liberty to choose the 
direction of spin change, since that can be allowed for by (3); we are also 
at liberty to choose M,, and we make our choice so as to give a minimum 
of calculation in the summation (1). We do this by setting J; equal to 
the larger of the initial and final spins and putting M;—J;. Then, 
beginning with M,=J;, the summation extends over increasing values of 
M ,— Mj, and selection rules operating on M;—M, reduce the summation 
to a few terms. When one of the spins involved is zero, the summation 
reduces to one term. 

For matrix elements such as 7’,,, S,;,, in which the operator combines 
large components of the initial nucleon four-component wave function 
with large components of the final one, relativistic wave functions are 
unnecessary, and we use two-component wave functions in the Pauli 
representation. If |+1) and |—1) are the eigenvectors of the operator 
corresponding to the z-component of intrinsic spin and belonging to the 
eigenvalues 4% and —43h respectively, then the individual particle wave 
functions are of the form 


b(n, l, m) [tueys b(n, l, m) [y) 
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“= 


where », J and m are the usual quantum numbers. However, in me 

i c i omponents 
matrix elements, for example A ;;, the operator combines large ¢ P : 
with small components, and so the wave functions used should be four- 
component solutions of Dirac’s equation. If xy, and x, are the large 
components of the solution, approximate values for the small components, 
xX, and x,, are given by :— 


il (lege oa 
X= 9M \\oe—* dy) X4° G2)? 


a Cua ) 
i= oat (= +t =) ee ats} 
where M is the ratio of the nucleon mass to the electron mass. 

If the nucleon energy is not much greater than its rest energy, x3 and 
x, are solutions of Schrédinger’s equation corresponding to the two 
eigenvalues of one component of the nucleon spin, and so we can extend 
the non-relativistic nucleon wave functions to four-component approxi- 
mate relativistic ones if we replace 


y(n, l,m) |+1) by :— and (n, l,m) |—1) by :— 
a Ov =i yf 2 eee 
Xa= OIL Bz X= oat ae —* dy) ¥ 
a | i Oe 
Xa= Of (F a4 a) X= — 90 Bz 
Xs= Xs=0 
X1=0 X= 


There are a few points to be considered in connection with the deter- 
mination of nuclear wave functions in terms of the individual particle 
wave functions. Firstly, we must decide which system of angular 
momentum coupling is likely to be the better approximation, the Russell— 
Saunders or RS form, or the j~j form. This will depend on the particular 
nuclei in which we are interested, the requirement being a plausible 
nuclear configuration whose wave function is an eigenfunction of angular 
momentum corresponding to a given eigenvalue. In dealing with shells 
of high / number which are about half-filled, the number of ways in which 
the particles can be accommodated is very large and there may be a few 
possible wave functions corresponding to a given pair of values J, M. 
These wave functions are much too unwieldy for calculations. In such 
cases it is simpler to use j—j coupling since this requires smaller shells, and 
we can make the simplifying assumption that particles can pair off to give 
zero contribution to the spin. 

Secondly, we are not interested in shells which are complete in the 
initial and final nuclei. The operator Q,, operating on a nucleon in a 
closed shell in the initial nucleus, gives zero if the nucleon is in a proton 
character state, and if it is in a neutron state it changes to a proton state, 
thereby becoming orthogonal to its opposite number in the final nucleus. 
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Finally, since the operator in the nuclear matrix elements operates on 
only one nucleon wave function at a time, we get contributions to the 
matrix elements only from those parts of the initial and final nuclear wave 
functions which are exactly the same except for, at most, one of the 
individual particle wave functions. In all other cases we get zero 
contribution because of the orthogonality properties of these wave 
functions. 

We consider now the matrix elements for the transition 36C143SA. 
The spin of this chlorine isotope has been measured by Townes and Aamodt 
(1949), who obtained the value 2, and 3A, being an even-even nucleus, 
is assumed to have zero spin. From the table given by Mayer (1949) a 
plausible configuration for *°Cl, assuming j~j coupling, would be :— 


Protons (1s)"(1p)°(1d5y5)°(1dg/9)*(2s)? 
Neutrons —_(1s)*(1/p)®(1dy/.)®(28)?(1dg))® 
For #®A the protons and neutrons both have the same arrangement, 


namely 
(1s)*(1p)°(1d5y2)°(2s)"(1dgy 5)”. 


Transitions from the chlorine arrangement to the argon arrangement 
conform to the parity requirements of second forbidden transitions. In 
accordance with the discussion above, we set J,;—2, M,=2, J;=0, M;=0, 
and so we must make up a wave function from the chlorine arrangement 
with J=2, M=2, and a wave function from the argon arrangement with 
J=0, M=0. For chlorine, we assume that the two protons in ldz/, 
states are paired together to give zero contribution to the nuclear spin, and 
we make the same assumption about the two protons in ldg), states in 
argon; thus the sub-wave-function of these two protons is the same in 
both nuclei and we need not consider them. For the remaining nucleons 
in unfilled shells in chlorine, the only wave function corresponding to 
M=2, taking into account the exclusion principle, is 

{é, 2, BN}. 2, a) NHK, 2, —$)N}{(2, 0, 4)P}. Sai (3) 
The WN and P indicate the character states of the nucleons and the numbers 
in brackets are the j, J, and m quantum numbers. This expression 


represents the Slater determinant giving a completely antisymmetrical 
wave function for the four nucleons :— 


(3,2,3)M(1) (, 2, 2)N(1) 
(3, 2,8) N(2) (¥, 2, )N(2) 


v 
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— 


(41-12 


For 36A, the two neutrons in 1d). states and the two protons in 2s states 
can combine in two ways to give M=0; these are 


A={(8, 2, $) N}(E 2, —§) HG, 0, YPHG, 0, —3)P} «  (44) 
B=({(3, 2, $) NG, 2, =) NG. 0, 4)PH{G, 0, —4$)P} . (4 b) 
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and the eigenfunction for J=0, M=0 is some combination of A and B. 
Examination of (3) and (4), however, shows that A gives no contribution. 
to the matrix elements since two of the individual particle wave functions 
are different in the initial and final nuclei. The whole contribution to the 
matrix elements comes from B, and since we are interested only in the 
ratios of matrix elements for the purposes of curve fitting, a multiplying 
constant for B is unimportant. Thus with (3) as the initial wave function 
and (4 6) as the final one, the matrix element of an operator O reduces to 


f{4, 0, —DP}OQ{(G, 2, Np da - =. ee} 
The angular and spin part of the wave function (4, 0, —$) is Yo9|—1) and 


for (3, 2, 3) it is 5-1/2(Y,, |+1)+2Y..|—1)). (Here the Y, ,, follow Rose’s 
definition (1937).) 

For the radial part of the wave functions we assume that the nucleons. 
move in an oscillator potential well, so that the particle wave function is a 
product of three functions, corresponding to the three co-ordinates, such 


Fo(e)—exp (—dle*). H, (When). |. 6) 


where & is a parameter such that k—1!/? is of the order of magnitude of the 
nuclear radius, and H,(€) are the Hermite polynomials. In the table 
giving her proposed scheme for nuclear shell structure, Mayer (1949) gives 
the spectroscopic term and the oscillator quantum number N corresponding 
to each nucleon state, and we combine products of the functions (6) to give 


x Cn gnymg Fn ,(*) Fay) Fn) . : : te (7): 


where n,+-7,+n,=WN ; the function (7) is the product of a radial function 
and a spherical harmonic corresponding to the given spectroscopic term. 
The radial function for the 2s states, with a corresponding spectroscopic 
term 2s and oscillator number 2, is 


4h3\ 1/4 
(=) exp (—thkr?) (2kr?—8). 


For the 1d states the oscillator quantum number is 2 and the spectroscopic 
term 3d, so that the radial part of the wave function is 


16\1/2 /k*\ v4 | 
' (35) (=) r exp (—dkr*). 


Thus the wave functions in (5) become 
1 1 4h3\ va : 
(5. Oe 5) = (=) exp (—$ghr?) (2kr?—3) Yoo [—15. 5. . | (8) 


5 ta 16\ 1/2 /k7\1/4 | ' 
ay (zs) @) ® exp (—$hr*) 5-1"(Yy1|+1)+2¥a9|—1)) 
(9) 
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and we are now in a position to work out for this particular disintegration 
those matrix elements which do not require extension of the wave functions 
to the four-component representation, namely f,, and 7';;. 


We get 


is y* Phe Ye 2H xy 
2 
if —l 0 0 0 1) Xo 
Wich 
so that the matrix elements R,, are 
Re= j 
—> 
1 } 0 
eee 
eee | V15k 
i) ade he 
and 
aes yf Ee We 2a ay 
oF 0 0 0 a 1 0 i 
Cy 0 0 0 —1 a O+} X VIBE’ 
ee es: 2 0 0 eee 23 


A check on these values is that the matrix elements S;,, should be zero, 
since they forbid the spin change 0<—>2. From these values we get 


Py= i 
‘| pe es | 0 r 
aon.) 0 0 


In order to evaluate the matrix elements A,, and A;,°, we must extend 
the wave functions in (8) and (9) to the four-component representation. 
Doing this we find that A,, and A,,’ are all zero for this transition. 

We evaluate now the nuclear matrix elements for the transition 


eT ch Ru. From Mayer’s table a plausible configuration for Tc is :—_ 


3 Protons in 1gg/.,__ 
6 Neutrons in 2d;,, (complete subshell), 


giving a spin of 9/2, in agreement with the experimental value (Kessler and 
Meggers 1950), and for *Ru 

4 Protons in 12/5, 

5 Neutrons in 2d;/., 


giving a spin of 5/2. Thus we get a spin change of 2 units and no oneeC 
in parity, conforming to second forbidden requirements. We set J,=9/2, 
M =9/2, and sum the transitions to J;=5/2, M,=5/2, 3/2. (We need go 
no further since 4M=4 is not allowed by any second forbidden matrix 
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elements.) The six neutrons in 2d,). states in *Tc form a complete 
subshell and so give no contribution to the nuclear spin, and, of the three 
protons in 1gy). states, two can be paired off to give zero spin. The 
remaining proton is therefore in a (§, 4, 3) state. For ®Ru the 4 protons 
are paired off to give zero spin ; for M,—5/2 all the 2d,,. neutron states are 
filled except (3, 2, —3), and for M,—3/2 they are all filled except (3, 2, —). 
We therefore evaluate the matrix elements of operator O for the two 
transitions, by considering 


on i 4, —)P} oa 2, —9)N} dr. 


For the 1g states, the oscillator quantum number is 4 and the spectroscopic 
term 5g, and so the wave function for the (2, 4, —$) state is 


128 \ V2 /f12\ 1 . 
(5) (=) r4 exp (—tkr?) Y,-4|—1). 


For the 2d states, the oscillator quantum number is also 4 and the spectro- 
scopic term is 4d, and so the wave function for the (3, 2, —5) state is 


8 \1/2 /&7\ 1/4 
(a5) (=) (2kr4— 7r*) exp (—3kr?) Vien | 215 


and for the (3, 2, —3) state is 


8 \1/2 /k7\ 1/4 
(was) (=) (2kr4— Tr?) exp (—4kr?)5 Y,-,|+1)—2Y,_ 1]—1)). 


7 


We are now in a position to evaluate the matrix elements R. 
For the transition to the M,—5/2 state, these are :— 


a 


ij? ij? 


k= 
il a 0 
. C 
7 —l 0) _ 
0) aC) 


alt ink: ae 
C being a constant arising from the normalization factors, 


Ly= 
47—l 0 
20 
=) () a, 
v) x bE? 
0) 0 0 
Siis= 1 


: 2C 
Sis. te xX rae all other S;,,=0. 
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For the transition to the M,=3/2 state, we get :— 


Sin= 1 

Sis= 7 } 

(oer x< Bh? all other S3,=0. 
122 

Soo2=—1t 


When we extend the wave functions to the four-component representation, 
we find that A,, and A,,’ are all zero for both transitions. 

The ratios of the matrix elements which we have found are summarized 
in table 2. Products containing R,,”, 7, are found from 


ik i=— T iz, 
and all other matrix elements and products which are not quoted are zero. 


ij? H) 


Table 2. Ratios of Matrix Elements 


2 | Ry? 2 |7;;\? 2 Sil? | 2 4° T;, 
0) 7) dik ij 
36(] 85364 1 9 0 3i 
9976 Ru 1 4 = 24 


It is worth adding here that, although these matrix elements have been 
evaluated with an explicit form of the radial wave function for the individual 
particles, the ratios obtained (i.e. the ratios of successive columns in table 
2) are substantially independent of this form. Certainly this is true for 
those matrix elements which combine large components of the wave 
function with large components, since the same integral 


{ Rk, Rr*dr 
0 


is involved in each case. It seems probable that the matrix elements 4A ,, 
would not be much affected by use of a different radial function, and 
indeed in cases where the final state of the disintegrating nucleon has the 
same radial function as the initial state, A,; is completely independent of 
the form of the function. The integrals involved are of the form 


if R?r2dr, which is 1 by normalization, 
0 


and im RR dr=—$ | Rr? dr=—3. 
0 0 


The detailed forms of the tensors 7';,, R ,;, etc., would be the same no matter 
what radial function were used, since all elements of a given tensor 
would be multiplied by the same factor. This is of particular importance 


in scalar products suchas 2 7';,*R,,, ete. 
0] 
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§ 5. COMPARISON WITH EXPERIMENT If 


Much of the arbitrariness involved in fitting the theoretical correction 
factors on to the experimental curve is eliminated now that we have 
estimates for the ratios of matrix elements. The problem is now to fit 


API+ AAgG+A7h 


on to an experimental curve e, where ,, Ag are constants determining 
the proportions in which two given interactions are mixed, if and h are the 
correction factors for the two interactions singly, and g is the cross 
correction factor; f, g and h are evaluated using the ratios of matrix 
elements. We have therefore to choose Aj, A, so that 


2 (APP AA ag HAg*h —e;)” 


is a minimum, the summation being over points on the curves. Minimiz- 
ation of this expression leads to two simultaneous third degree equations 
for A, and A, which would be rather difficult to solve. However, it was 
pointed out in §3 that approximate linear relationships frequently exist 
among the coefficients of the matrix elements, and we make use of this 
property here by looking for a linear relationship between f, g and h, i.e. 
we express / in terms of f and g, 


h=af+ Bg. 


The problem of curve fitting then reduces to the fitting of two curves, 
i.e. we now make 


» (af +69 ;—e;)? . 


a minimum, and the values of a and b, along with « and £, give us A, and Ag. 
In all the cases considered here, such linear relationships have been found, 
almost all of them accurate to within 1%. 

The beta spectrum of ?°Cl has been measured by Wu and Feldman 
(1949), using a source of thickness 0-3 mg/cm?, and they found that when 
the Fermi plot is corrected by (S,;; shape) the result is a straight line. 
However, the measured spin of ®°Cl is 2 and so this cannot be the explan- 
ation of its spectrum. Longmire, Wu and Townes (1949) tried to fit the 
spectrum with a second forbidden correction factor using arbitrary values 
for those matrix elements which allow a spin change of 2-+0. They found 
this impossible to do, and went on to try the effect of a combination of 
interactions. According to them it is possible to fit the spectrum with the 
combinations S7', V7’ or VA, but. not with SA. Since then, Wu 
and Feldman (1951) have revised the attempt to fit the spectrum with a 
pure tensor interaction, and found that it is possible to do so. Similar 
results have recently been obtained by Fulbright and Milton (1951), 
whose measurement of the spectrum agrees with that of Wu and Feldman, 
except that they found fewer electrons at low energies. 
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The measurements of Wu and Feldman have been used here for 
comparison of theory with experiment. The process of curve fitting has 
been carried out for the four combinations of interactions as before and 
also for the tensor interaction by itself. (In this case there is only one 
curve for the tensor interaction.) None of the interactions gives a really 
close fit. Certainly the tensor interaction, ST and VA can be ruled out. 
The curves obtained from SA and V7 are very similar and are fairly close 
to the experimental curve. They are shown in fig. 4(a). However, the 
coefficients obtained in fitting the correction factor for S and A give 
imaginary values for the ratio Ag: A4, and since we assumed these constants 
to be real, we cannot accept this fit. For V and 7 we obtain 


The corrected Fermi plot, obtained by dividing the original by the square 
root of the theoretical curve in fig. 4 (a), is shown in fig. 4(b). It deviates 
from a straight line at low energies, but above this it is almost straight. 


Fig. 4 


N 
Fe@(Wy-W)? 


1:0 [Sait 2-O 2:5 a |:0 pS as 2-0 2-5 


(a) (ON 
360], Best fit for V7’, using calculated *6Cl. Corrected Fermi plot using 
ratios for matrix elements. correction factor shown in fig. 4(a). 


It is interesting to note that this deviation is in the direction which would 
result from source thickness distortion, and that it would be reduced using 
Fulbright and Milton’s results. oo . 

This process has been applied also to the transition ™Tc%S*°Ru, 
using the ratios of matrix elements in table 2. Both SA and V7 give good 
fits, as shown in fig. 5(a), but the others are very poor. The ratios of the 
constants obtained are 


AS 1.5 or 0:66, 
Aa 
Ay 1.6 ot 0-63, 
Ap 
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and if we take this in conjunction with the results for °°Cl and the results 
of §3 we conclude that the interaction is approximately two parts of vector 
to one of tensor. When one considers the possibility of experimental 
errors, the uncertainty in the nuclear radius, and the crudeness of the 
nuclear model used for the evaluation of the matrix elements, the agree- 
ment between these results is very satisfactory. The corrected Fermi plot 
of Tc for the V and 7 combination is shown in fig. 5 (6); it is a straight 
line for most of the energy range, deviating from it only at low energies, 
again in the way which would be caused by distortion due to source 


thickness. 
Fig. 5 


1-0 1-2 w V4 1-6 


°°Tc. Best fit for VZ' or SA, using Tc. Corrected Fermi plot using 
calculated ratios for matrix elements. correction factor shown in fig. 5(a). 


§ 6. CONCLUSION 


In recent months there has been a number of papers published in 
connection with a combination of the five interactions in beta-decay. 
Moszkowski (1951) has applied matrix elements calculated by Wigner 
from group-theoretical considerations to the difference in ft values of 3H 
and ®He, and also to the branching ratio in the decay of 7Be, and obtains 
the result that the interaction is approximately 1 part of Fermi type 
(S or V) to 2 parts of Gamow-Teller type (7' or A). The cross correction 
factors for an arbitrary degree of forbiddenness have recently been 
published by Pursey (1951 b). He applies his results for first forbidden 
transitions to several nuclei whose spectra have the allowed shape but 
whose ft values are high, and concludes that it is possible to explain these 
as first forbidden transitions if the interaction is mainly 7’ with a small 
admixture of S or V, or mainly V with an admixture of A. 

Biedenharn and Rose (1951) have shown that the assumption of real 
values for the A in the invariant is no restriction on the theory. 

From a detailed consideration of ft values for allowed transitions, Horie 
and Umezawa (1951) dedude that the interaction is mainly Fermi type 
(S or V) and that the Gamow-Teller type (7’ or A) is subsidiary. Finally, 
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Tolhoek and de Groot (1951) have proposed two principles of symmetry 
in consequence of which the coefficients \ in the invariant must be real 
and the interaction can be combinations of 8, A and P only, or a combin- 
ation of V and 7 only. 


In conclusion, then, our result that the interaction is two parts vector 
to one of tensor disagrees with the considerations of Moszkowski and of 
Pursey, but agrees with those of Horie and Umezawa, and it conforms to 
the requirements of the symmetry principles of Tolhoek and de Groot. 
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SUMMARY 


The energy distributions of the alpha-particles emitted from small 
nuclear disintegrations caused by the cosmic radiation in photographic 
emulsion have been found. These are analysed into the parts due to 
the light and the heavy elements of the emulsion, and the variation of 
nuclear temperature and effective potential barrier height with star size 
found for each. 


§ 1. INTRODUCTION 


SEVERAL investigators (Harding, Lattimore and Perkins 1949, Perkins 
1950, Bernardini, Cortini and Manfredini 1950, Cortini 1950) have 
measured the energy distribution of the alpha-particles emitted from 
nuclear disintegrations in photographic emulsions exposed to the cosmic 
radiation. These results give the mean energy distribution averaged 
over stars of several numbers of tracks, or sizes. This tends to obscure 
any variation of the energy distribution with energy. 

In this paper, the alpha-particle energy distribution is given for stars 
of several sizes. The contributions of the disintegrations in the light and 
heavy elements of the emulsion are evaluated separately, and the variation 
of the nuclear temperature and potential barrier height with star size 
evaluated for each. 


§ 2. EXPERIMENTAL PROCEDURE 


The observations were made on disintegrations in C2 emulsion exposed 
horizontally to the cosmic radiation on the Jungfraujoch and processed 
in the usual manner. C2 emulsion was chosen rather than the more 
sensitive G5 so that the tracks of alpha-particles could be readily 
distinguished from those of protons. 

The stars were classified according to their numbers of tracks and, for 
those of less than seven tracks, according as they had a recoil or not. 
A recoil is a short, thick track due to the residual nucleus. Tracks less 
than four microns in length were counted as recoils. 

One thousand stars were examined, and the ranges of all the 
alpha-particles emitted from them which ended in the emulsion were 
recorded, The energy of each alpha-particle was found using the 
SE eee 
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range-energy relations of Lattes, Fowler and Cuer (1947), extrapolated 
when necessary to higher energies. The usual geometrical corrections 
to allow for those particles leaving the emulsion were applied throughout. 


§3. EXPERIMENTAL RESULTS 


The energy distributions of the alpha-particles from three-, four-, and 
five- and six-track stars, with and without recoil fragments, are shown 
in figs. 1 to 6. 


Fig. 1 


Histogram—Observed 
Distribution 


Fitted evaporation curves : 

Light elements 
— — — Heavy elements 
—— Total 


Corrected number of alpha-particles per 2 Mey. 


0 5 10 15 20 25 30 35 40. 
Energy in Mev. 


Energy distribution of alpha-particles from three-track stars with recoil: 
fragments, together with fitted evaporation curves. 


Comparison of the energy distributions for stars of the same size but 
with and without recoil fragments enables the contributions of the stars. 
in light and heavy elements to be evaluated separately. 

Harding (1949 b) has shown that most of the stars with recoil fragments. 
are due to the disintegration of the heavy elements silver and bromine, 
while most of the stars without recoil fragments are formed in the light 
elements carbon, nitrogen and oxygen. ‘There are therefore more 
alpha-particles from heavy elements among those from stars with recoils, 
and less from those without recoils, than from all stars taken together. 
The following procedure was used to separate the total energy distribution 
into its component parts due to disintegrations in the light and heavy 
elements. It is based on the assumptions that there are no low energy 
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(<8 ev for three-track stars ; <6 Mev for four) alpha-particles emitted 
from heavy elements and no high-energy (>15 Mev) alpha-particles 
emitted from the light elements. 


Fig. 2 


Histogram—Observed 
Distribution 


Fitted evaporation curves : 
Light elements 
— — — Heavy elements 


Corrected number of alpha-particles per Mev. 


0 5 10 15 20 25 30 35 40 
Energy in Mey. 


Energy distribution of alpha-particles from three-track stars without recoil 
fragments, together with fitted evaporation curves. 


Fig. 3 
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Histogram—Observed 
Distribution 
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Fitted evaporation curves: 

Light elements 
— — — Heavy elements 
———- Total 


Corrected number of alpha-particles 
per 2 Mev. 
— 
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0 5 


10 15 20 25 30 35 40 
Energy in Mev. 


Energy distribution of alpha-particles from four-track stars with recoil 
fragments, together with fitted evaporation curves. 
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Fig. 4 


Histogram—Observed 
Distribution 


Fitted evaporation curves: 

Light elements 
— — — Heavy elements 
—-—— Total 


Corrected number of alpha-particles per Mev. 


0 5 10 15 20 25 30 35 40 
Energy in Mev. 


Energy distribution of alpha-particles from four-track stars without recoil 
fragments, together with fitted evaporation curves. 


Fig. 5 
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Distribution 
——-— Fitted evapora- 

tion curve 
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Corrected number of alpha-particles per Mev. 


0 5 10 15 20 25 30 35 40 
Energy in Mev. 


Energy distribution of alpha-particles from five- and six-track stars with recoil 
fragments, together with fitted evaporation curve. 
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The observed energy distribution for stars without recoils was in each 
case multiplied by the factor which made the number of particles of low 
energy (as defined above) equal to the number of particles in the same 
energy region of the corresponding distribution of particles from stars 
with recoils. The energy distribution of the alpha-particles from the 
heavy elements alone was then obtained by subtracting these two 
distributions. This was fitted by the best evaporation curve of the form 

P(E) dE= eal exp (- me | dE 
where 7’, the mean nuclear temperature, and V, the mean potential 
barrier height for the de-excitation process, are adjustable constants. 


Fig. 6 


Histogram—Observed 
Distribution 


Fitted evaporation curves : 
Light elements 
— — — Heavy elements 


Total 


Corrected number of alpha-particles per Mev. 


0 5 10 15 20 25 30 35 40 
Energy in Mev. 
Energy distribution of alpha-particles from five- and six-track stars without 
recoil fragments, together with fitted evaporation curves. 


The energy distribution for five- and six-track stars in heavy elements 
was obtained by fitting an evaporation curve directly to the distribution 
for stars with recoils, as it is most unlikely that any of these could have 
been formed in the light elements. 

Owing to the small probability of an energetic alpha-particle ending 
in the emulsion, the correction factors of those that do are rather high. 
This makes the final values of 7' and V rather sensitive to the number 
of such particles in the distribution concerned. In order to minimize 
this effect, the distribution for energies above 20 Mev was added for all 
stars of three to six tracks, both with and without recoil fragments, and 
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divided among the individual distributions in proportion to the number 
of particles from heavy elements as determined from a preliminary 
calculation without high-energy averaging. This procedure does not 
invalidate the results, since Bernardini, Cortini and Manfredini showed 
that the distribution of the high-energy alpha-particles varies little, if at 
all, with energy. The distributions shown in the figures are as observed, 
but the final values of 7 and V were calculated using this high-energy 
averaging. 

These energy distributions of alpha-particles from stars in heavy 
elements were then subtracted from the corresponding energy distributions 
for stars without recoils, so that the number of particles of more than 


Fig. 7 


< Mean nuclear temperature 
O Mean potential barrier 

~ height 
——— Le Couteur (1950) 


° 


T or V in Mey. 


Star size. 


Variation of mean nuclear temperature and mean potential barrier height 
with star size for disintegrations of heavy elements. 


15 Mev was reduced to zero. The energy distributions of the alpha-particles 
from the light elements obtained in this way were also fitted by evaporation 
curves, and the best values of 7’ and V determined. 
_ The evaporation curves found in this way are included in figs. 1 to 6 
in each case, and the corresponding values of 7’ and V are plotted in 
figs. 7 and 8 as a function of star size for the heavy and light elements 
respectively. 

The energy distribution of alpha-particles from one-track stars formed 
by 7--mesons given by Menon, Muirhead and Rochat (1950) has been 
analysed in a similar way, and the resulting values of 7 and V are also 


plotted in figs. 7 and 8. 
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The alpha-particles from the disintegrations of light elements may 
also be found using the sandwich emulsion technique developed by 
Harding (1949 a). A composite emulsion consisting of alternate layers 
of pure gelatin and emulsion was exposed and processed. It was 
examined for disintegrations in the gelatin layer, which must occur 
in carbon, nitrogen or oxygen. Twenty-six such disintegrations, of 
mean size 3-65, were measured, and the energy distribution of the 
alpha-particles found. The values of 7 and V appropriate to this 
distribution are also plotted in fig. 8. 

The values of the potential barrier height for the unexcited nuclei are 
interpolated from results given by Bethe (1937). 

As these energy distributions were fitted by a formula having two 
adjustable constants, the values of Z’ and V are subject to some 
uncertainty, which is estimated to be of the order of 30%. 


Fig. 8 


x Mean nuclear temperature 
O Mean potential barrier 
height 


T or V in Mev. 


Star size. 


Variation of mean nuclear temperature and mean potential barrier height 
with star size for disintegrations of light elements. 


This method of analysing the total energy distribution into its 
components due to disintegrations in the light and heavy elements 
rests on the assumption that there are no alpha-particles of less than 
8 and 6 Mev emitted from three- and four-track disintegrations in heavy 
elements. This is, however, rather less stringent than that of Perkins 
(1949), who assumed that all alpha-particles of less than 10 Mev are 
emitted from light nuclei. 

It was also assumed that there are no high-energy alpha-particles 
emitted from light elements. The results from the sandwich plates 
showed that such particles are sometimes, though rarely, emitted. 
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As stars of the same size are produced by particles of a range of 
energies, the results for each star size correspond to the mean of a wide 
distribution of excitation energies. This mean excitation energy, however, 
increases regularly with star size. 

The energy distributions of the alpha-particles emitted from dis- 
integrations of more than six tracks were also measured, and found to 
be in good agreement with the results of Harding, Lattimore and Perkins 
(1949) and of Bernardini, Cortini and Manfredini (1950). 


§ 4. Discussion 


The results given in figs. 7 and 8 show that, for the disintegrations of 
heavy nuclei, the nuclear temperature rises and the mean effective 
potential barrier falls with increasing star size which corresponds, on 
the average, to increasing excitation energy. These results are in accord 
with expectation, since the more highly excited a nucleus is, the higher 
its temperature and the more its effective potential barrier is reduced 
by its violent oscillations. The variation of nuclear temperature with 
star size is in good agreement with the calculations of Le Couteur (1950), 
which, however, refer to the initial temperature. But since the alpha- 
particles are probably emitted at the start of the disintegration, this is 
not much different from the mean temperature over the period of 
alpha-particle emission determined here. 

The results for disintegrations in light nuclei show similar trends for 
small star sizes, but both the temperature and the barrier height appear 
to level off to constant values for the more energetic disintegrations. 
This is what would be expected if there is a limit to the excitation energy 
of a nucleus, and hence also to its temperature. The de-excitation of 
a nucleus excited to this maximum temperature results in its complete 
disintegration. The potential barrier is very low during the latter part 
of this process and so its mean value over the whole de-excitation is but 
a small fraction of its value for the nucleus in its ground state. This 
maximum excitation is only attained for the heavy nuclei at much 
higher energies than those investigated here. 
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XCI. Some Viscous Fluid Flow Problems 
I: Jet Emerging from a Hole in a Plane Wall 


By H. B. SQuirE 
Imperial College, London, 8.W.7* 
[Received May 30, 1952] 


SUMMARY 


An exact solution of the equations of viscous fluid flow for axially- 
symmetric motion is studied for a special case : this is shown to give the 
flow for a jet emerging from a hole in an infinite plane wall. Some 
examples of the calculated streamlines are given. 


§1. [INTRODUCTION 


It was shown in a previous paper (Squire 1951) that some exact solutions 
of the equations of viscous fluid flow for axially-symmetric flow could 
be obtained by assuming that the stream function ¢% is of the form 


we=vrf(e) > ne ne 


where 7, 0, ¢ are spherical polar coordinates with @ measured from the 
axis of symmetry (fig. 1), ~=cos @, and v is the kinematic viscosity. The 
component velocities are (u,v, 0) measured in the directions (r, 6, d) 
respectively, and 


Leow ape 
u= and 06 rt WH) 
Lo LO er 
| ee ee (2) 


~ rsing dr ra/(1—p2) 
These velocities are inversely proportional to the distance from the 
origin. It was also shown by Squire (1951) that the function / satisfies 
the equation 


P= 4uft-2(L— yt) pS) ee (B} 


where (4) is a quadratic function of whose three coefficients are 
arbitrary, being constants of integration of the fourth order differential 
equation satisfied by f. The solution for the special case S(.)= 0 was 
studied and shown to give the flow for a jet in an infinite expanse of 
fluid. 

A further example of the solutions of (3) is derived in the present note. 
gn ne 


* Communicated by the Author. 
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- §2, TRANSFORMATION OF THE EQuaTIon FOR THE STREAM FUNCTION 


Equation (3) can be transformed into a linear equation by a standard 
method such as that given by Forsyth (1929). Putting 


[} 


fa—21—w)E SAMAR. Steet AAA) 
and substituting in (3) gives 

»___2(u) 

g etait see iy ay esl oe oe (5) 


This equation is a form of the hypergeometric equation and known 
methods of analysis can therefore be applied. When two independent 
solutions g, and g, of (5) have been obtained the function f is given by 


Forsyth (1929) 
. aera 
= — 2(1—p?) | ——— om see bs erm ( 
f=—21—p*) | Se () 
where ¢ is a constant of integration. 


§3. SOLUTION OF EQuATION (5) IN A SPECIAL CASE 
If we take 4) to be of the form 
X(w)= — (40?+1)(1—p)? 
where b is an undetermined constant, (5) becomes 
Gq == 
The solutions of this equation are 
g=(1-+p)” 


(MEER. Ol 37) 


where 
n=$-+(—1)'0. 
It will appear later that b is real for the jet problem under consideration 
and then two (real) independent solutions of (7) are 
g,=(1+p)'? cos M, 
and 
go= (1+ pH? sin M, 
where 
M=b log (1-4). 
From (6) we then obtain 


sin M—c cos 7m) |: (8) 


ie | 1425 (a M-+csin M 


$4, Jet Emprcine From A HOLE IN A PLANE WALL 


We shall now show that the above solution can be applied to calculate 
the flow due to a jet emerging normally to a plane wall from a small 
hole. The wall is taken to be the plane p=0 (9=7/2). The boundary 
condition of zero velocity normal to the wall is satisfied but not the 
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condition of zero tangential velocity. The wall is a streamline and hence 
f=0 on p=0; also on the wall M=0 and hence from (8) 

1+ 2cb=0. 
Thus (1+.4b2)(1—p) (2) 

f= 2b cot [b log (1+y)J]—1° 
A numerical study of (9) shows that jet flow is obtained for positive 
values of 6 up to a value for which the denominator is zero on the axis 
of the jet »=1 (0=0). This occurs when 
2b cot (b log 2)=1 

for which the significant root is 

b=1-8937. 


Fig. 1 


ERAS aN aa 


PCS 


Jet emerging from hole in wall. 
b=1°85. F/pv?~1:2x 10°. 


The interesting cases of narrow (high speed) jets are obtained by taking 
values of the constant b just below this limit. 
Some examples are given in figs. 1-3 which show the streamlines for 
the cases 
b=1-85, 1-88, 1-8937. 
Figure 3 shows the inflow in the limiting case of a very high speed jet. 


§5. Discussion 


The flows calculated above and depicted in figs. 1-3 can be described 
in a general way as the fields of flow due to a jet emerging from an orifice © 
in a wall, for a range of values of the viscosity or of the jet momentum. 
There is, however, zero mass flow across the wall u.=0 arid the flows are 
more accurately described as the consequence of the application of a 
force F at the origin acting in the direction of the resulting jet (see 
Squire 1951). In general the determination of the magnitude of this 
force is rather complicated, but in the case of the narrow jets the force 
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1s approximately equal to the rate of change of momentum of the fluid 
crossing a plane parallel to the wall. The values of F/pv? determined in 
this way are given in figs. 1-3. 

It is interesting to note that the inflow into this jet is a viscous flow 
and is not exactly the inflow towards a uniform line sink situated on the 


eee” which would at first sight be expected: this is clearly shown 
in fig. 3. 


Jet emerging from hole in wall. 
b=1°88. F/pv? ~3°8 x 103. 


Fig. 3 


Limiting case of inflow into high speed jet. 
b=1:8937. F—o. 


§6. FuRTHER WORK 
The flow due to a line sink normal to a wall with zero velocity at the 
wall is included among the solutions of (3) and (5); this problem is now 
being investigated. 
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XCIL. Heat Conductivity of Superconductive Lead below 1° K 


By J. L. Otsen* and C. A. RENTON 
Clarendon Laboratory, Oxford+ 


[Received June 21, 1952] 


Ix view of the particular interest which the heat conductivity of 
superconductors has at very low temperatures the measurements by 
Mendelssohn and Rosenberg (1952) on a lead single crystal have now been 
extended into the range below 1° K. Measurements on tin and tantalum 
at very low temperatures have been reported by Daunt and Heer (1949) 
and Goodman (1951) but no directly observed data have so far been given. 
In both these cases the temperature was measured at the ends of the 
specimen by pills of compressed paramagnetic salt of high heat capacity. 
In order to avoid the uncertainties in interpretation involved in this 
method we have employed the standard procedure of measuring the 
temperature at two intermediate places along the specimen by means 
of carbon thermometers. Heat was supplied electrically at one end of 
the specimen while the other was thermally connected to a pill of 
paramagnetic salt which merely served as a heat sink. 

The heating wattage was of the order of 10u watts. This may be 
compared with the heat dissipated in the thermometers which lay 
between 0-1 watt at 1° and approximately 0-02 watt at the lowest 
temperatures. The heat influx through the leads amounted to somewhat 
less than 0-1, watt. The small mass of the carbon thermometers and 
the low specific heat of the specimen made the attainment of temperature 
equilibrium very rapid, and permitted reliable calibration of the 
thermometers. This also meant that the thermal conductivity of the 
specimen during a magnetic cycle could be measured without the 
temperature of the pill changing by more than a few hundredths of a 
degree between the beginning and the end of the cycle. (There is of course 
always some magnetic heating when a field is applied, and cooling when 
it is removed.) 

RESULTS 


The thermal conductivity in the superconducting state is shown 
plotted against the temperature in fig. 1. The normal conductivity was 
too large to allow measurements on the same specimen, but the results 
of Mendelssohn and Rosenberg (1952) on this sample indicate that it 
would be of the order of 6 watt/degree cm at 1° K. 

Figure 2 shows a logarithmic plot of the results yielding a rough 
proportionality of the heat conduction to 73 between 0-4° K and 0-9° K. 
Above this temperature the rise becomes more rapid. If, as might be 
expected in the present case (Olsen 1952), the electronic scattering of } 
the lattice vibrations is very small and if the heat conductivity is purely 


* Now at the Eidgenéssische Technische Hochschule, Ziirich. 
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due to the lattice then according to Casimir’s (1938) calculation it should 
be equal to AZ%. However, it should be noted that the theoretical 
value of the constant A is five times higher than the observed one in 
this case. 

In fig. 3 the heat resistance of the specimen in 3 magnetic cycles at 
temperatures of 0-43° K, 0-59° K and 1-5° K is shown. In all three cases 
there occurs the rise of heat resistance in the intermediate state observed 
earlier by Mendelssohn and Olsen (1950) and also by Webber and Spohr 
(1951). In the two cycles below 1° a correction has been made for 
the slight temperature change caused by the application of a magnetic 
field to the pill. 


Fig. 3 


deg em/watt 


Heat Resistance. 


vee 
Magnetic Field. 
Variation of heat resistivity with magnetic field at 0-43°K, 059° K and 


135° KO O Field increasing, @- — —@ Field decreasing. In all 
cases the heat resistance in the normal state (H >H,) was too small to be 
measured. 
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XCIII. The Growth of Crystals of the Hexagonal Metals from their Vapours 


By A. J. Forty 
H. H. Wills Physical Laboratory, University of Bristol* 


[Received June 5, 1952] 


ABSTRACT 


The techniques for preparation of well-developed growth surfaces of 
crystals of the hexagonal metals are described. A subsequent examination 
of the surfaces under the phase contrast microscope reveals growth steps 
which are probably monomolecular. Spiral growth hills on basal surfaces 
of magnesium and cadmium crystals indicate that these grow by the 
screw dislocation mechanism. The growth patterns on basal surfaces 
of zinc can only be explained by a spreading of (0001) monolayers, which 
are catalyzed by a foreign nucleus, probably zinc oxide. Similar nuclei 
incorporated in zinc crystals during growth are revealed by evaporation 
from the habit faces. The possibilities of performing plastic deformation 
experiments with these crystals are discussed. 


§1. GENERAL INTRODUCTION 


THE concepts of the Edge Dislocation (Taylor 1934) and the Screw 
Dislocation (Burgers 1939) were originally proposed, on purely theoretical 
grounds, to account for the mechanical properties of real metal crystals. 
There is already considerable experimental evidence in support of the 
application of the concept of the screw dislocation to account for the 
growth of real crystals from low degrees of supersaturation (Frank 1949, 
1952). The studies of growth surfaces, to be described in a later section, 
verify that with the possible exception of zinc, the hexagonal metals 
grow by the screw dislocation mechanism. This provides the first direct 
observation of the existence of screw dislocations in metal crystals. 

Of the possible methods for preparing single crystals of these metals, 
growth from their vapours in vacuo or in an inert atmosphere is the most 
satisfactory for this investigation. This provides crystallographic 
surfaces, suitable for a sensitive microscopic examination, which are 
relatively free from chemical overgrowths. 

A high ratio of ambient vapour pressure to equilibrium vapour pressure 
is necessary to nucleate stable crystals. The vapour pressures at the 
melting points, that is the triple point pressures, are listed for the 
common metals in table 1.¢ To nucleate crystals on the walls of a 
containing vessel it is generally necessary, initially, to heat a stock of 
metal above its melting point. The crystals are then developed at a 


* Communicated by the Author. ioe. 
+ Nore: These values are obtained from data given in the International 
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lower supersaturation by decreasing the temperature of the stock to 
below its melting point. For the hexagonal metals, whose triple point 
pressures lie between 10-1 and 10 millimetres of mercury, crystal 
platelets can be established on cooler parts of the stock without the 
initial overheating. The subsequent rate of growth of the habit faces, » 
and therefore the lateral separation of growth steps on these, depends on 
the degree of supersaturation close to them. This is governed by the 
shape of the vapour pressure curve in the region of the melting point and 


Table 1. Vapour Pressures of the Common Metals at their Melting Points 
(p in mm Hg) 


Melting Point 
Tu SAK 


1515 
692 
594 

2895 

1800 

1123 

1083 
932 

1233 
904. 

1728 

3660 

1336 
234 

1357 
370 
335 
600 


505 


the difference between the temperatures of the stock of subliming metal 
and the growing crystal. A supersaturation of about 0-001 is required 
for the growth steps to be easily resolved under the optical microscope ; 
that is, p/p, must be 1-001, where p,, p, are the equilibrium vapour 
pressures of stock and growing crystal respectively. With magnesium, 
for example, if the stock is at a temperature of 620° c the growing crystal 
must be only 1/30 degree lower. It might be considered that this 
necessitates careful control of temperatures, but platelets several 
millimetres wide with surfaces true to the crystallographic plane to within 
about 10 seconds of are have been grown with comparatively crude 
methods of heating. A very careful measurement and control of actual 


temperatures of the stock and growing crystals would be extremely 
difficult. 
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The essential feature of the techniques described below is the use of a 
single source of heat for both stock and crystals ; good crystals grow where 
there is a suitably small difference of temperatures between them and the 
stock. The constancy of temperature to be maintained by the heater 
depends largely on the shape of the vapour ‘pressure curve close to the 
working temperature. For magnesium, a tolerance of --10 degrees, giving 
only 5% deviation in the degree of supersaturation, is permissible. 


§2. PREPARATION OF THE CRYSTALS 


Flat hexagonal platelets of magnesium have been grown by subliming 
pure magnesium turnings (99-99%, Mg) at 620° in a slow stream of argon. 
The turnings are heated in a nickel foil liner in an electrically-heated 
nickel tube. Crystals, several millimetres wide, grow on the stock after 
about 15 hours. Those crystals growing on the outside of the stock 
usually have higher quality surfaces than those inside. 

Zinc has a triple point pressure higher than that of magnesium. Large 
flat hexagonal platelets can be grown by heating small pieces of zinc rod 
at about 400° c for five hours in a sealed quartz tube, initially evacuated 
to about 10°? mm-mercury. The larger crystals grow on a substrate of 
zinc deposited on the wall of the quartz tube when heating was 
commenced. If the metal is heated in a slow stream of argon at the same 
temperature, similar crystals can be grown in about 10 hours. 

It has been found that large basal faces of zinc can be prepared by 
heating the metal (99-9°% zinc) in argon or im vacuo in a quartz tube to 
well above its melting point. The zinc evaporates rapidly and condenses 
in large spherical droplets on the upper, cooler parts of the tube. When the 
heater is removed, the droplets cool slowly and large facets are developed 
on the free surfaces. These facets have a crystallographic quality 
comparable with that of the platelets grown directly from the vapour. 
It is thought that the perfection of the facets can be attributed to the 
freezing of a melt with a large free surface, and subsequent evaporation 
during the slow cooling of the solid. 

Cadmium has a lower triple point pressure than zinc and consequently 
it is rather more difficult to grow large crystals. Flat hexagonal platelets, 
suitable for microscopic examination, can be grown by heating cadmium 
rod or turnings in a sealed, evacuated quartz tube at 300° c for 24 hours. 
As in the case of zine, it is convenient, for easy nucleation of crystals, to 
deposit a substrate of cadmium on the walls of the tube by overheating 
initially. Basal surfaces of cadmium can also be developed by the 
droplet technique described for zinc. These are usually covered with 
cadmium dendrites and are unsuitable for the present investigation. 


§ 3. EXAMINATION OF GROWTH SURFACES 
The flat platelets, with the basal (0001) surfaces well developed, are 
very suitable for examination under the optical microscope. Usually, 
the surfaces depart from the true crystallographic (0001) plane by only 
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a few seconds of arc, so that a sensitive examination under reflection, 
phase constrast illumination is possible. The surfaces have a high 
reflectivity and the silvering technique normally adopted for the 
observation of monolayers on other crystals is unnecessary. The crystals 
are seldom bigger than 2 or 3 millimetres wide and are mounted on 
plasticine on a stage whose tilt can be adjusted so that their surfaces are 
normal to the axis of the microscope. Very small steps, less than 
20 A high, can be detected only when this condition is fulfilled. Under 
favourable conditions of illumination, depending on the absence of light 
scattered from nearby coarse growth structures or other crystals, growth 
layers which are almost certainly monolayers have been detected. 


§4. GrowTH PaTrERNS ON MAGNESIUM CRYSTALS 


The growth steps, illustrated in the following micrographs are barely 
visible through the visual eyepiece of the microscope. By employing 
suitable photographic techniques, the forms of growth hills can usually 
be seen on a micrograph. The low visibility of the growth steps suggests 
that their height is less than 20 A. (The phase contrast microscope used 
has been calibrated by measurement of visibilities of steps on silvered 
cleavage surfaces of mica, whose heights were measured by multiple-beam 
interferometry—Forty 1952.) Other considerations, such as the absence 
of a hollow core for the screw dislocations observed and the similar 
visibilities of steps in different growth hills on different crystals, indicate 
that the growth layers are monomolecular. 

Often, after the magnesium crystals have been mounted on the 
plasticine for some time, the surfaces become coarsely speckled. The 
speckle is frequently concentrated along the edges of the monolayers 
giving them a visibility of practically unity under phase contrast 
illumination. Figure 1 (Plate LXI) is a phase contrast micrograph 
of a small crystal with the growth steps ‘decorated’ in this striking 
manner; the steps appear bright on a darker speckled background. 
This particular crystal had been in contact with air and plasticine for 
24 hours. It is thought that the decoration is due to the creep of a film 
of reactive material, possibly fatty acid, from the plasticine. This might 
combine with magnesium oxide, which is expected to be concentrated 
along the edges of growth layers, to form a ridge of material at the steps. 
The phase contrast microscope is particularly sensitive to this form of 
surface singularity. The observations made on this phenomenon suggest 
that the exposures of the surface to both air and plasticine are critical. 
So far, all experiments to produce the ‘ decoration’ under controlled 
conditions by spreading films of stearic acid, oleic acid and liquid paraffin 
have been unsuccessful. This is not the first time that a ‘ decoration ’ 
of growth steps has aided their observation. Micro-etch pits concentrated 
along the edges of monolayers on natural beryl (Griffin 1951) and ridges 
of crystalline overgrowth, probably quartz, concentrated along growth 
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steps on silicon carbide crystals (Vand 1951) have greatly improved their 
visibility. 

Monolayers have been detected on about 50% of the crystals examined. 
It is a characteristic feature of most of the magnesium crystals that there 
is only one growth hill on a basal surface. In almost every case where the 
centre of a growth hill can be resolved, the step-line is a spiral or several 
spirals, centred on screw dislocations with Burgers vectors having a 
component normal to the basal plane. This is the first direct observation 
of screw dislocations in metal crystals; with the exception of a rather 
abnormal example in gold (Amelinckx 1952). 

Crystals growing on single screw dislocations are illustrated in figs. 1, 2, 3 
(Plate LXI). The step-lines are usually curvilinear but they tend 
sometimes towards a trigonal shape in accordance with the symmetry of 
the (0001) surface. 

Some growth hills consist of several parallel spiral ramps ; the step-lines 
terminate on a row or cluster of similar screw dislocations. Figure 4 
(Plate LXII) shows a crystal growing from a cluster of co-operating 
screw dislocations. Figure 5 (Plate LXII) illustrates a growth hill 
with eleven parallel spiral step-lines centred on a row of eleven equal 
co-operating screw dislocations. These rows or clusters of equal screw 
dislocations could result from the disintegration of a large screw 
dislocation in the platelet. There is, however, very little conclusive 
evidence for the formation of large screw dislocation groups in the 
hexagonal metals. 

Crystals growing from two or more screw dislocations of opposite 
hand have not been observed. There are a few examples of growth hills 
with step-lines in closed loops but there is no terrace at the centre ; 
the growth layers are probably nucleated by a speck of magnesium oxide 
or other foreign material in these cases. Figure 6 (Plate LXIT), 
which illustrates such a growth hill, is also of interest because it illustrates 
the domination of several screw dislocations by the growth layers from 
the foreign nucleus. 

Sometimes growth layers are nucleated from the edges or corners of 
the crystals. If the degree of supersaturation is high enough and the 
crystal contains no screw dislocations or foreign surface nuclei, then the 
corners of the surface are favoured sites for nucleation of layers. In some 
of the examples of growth from the corners, and probably in all cases of 
growth from the edges of the surface, contact with an adjacent crystal 
may be responsible for the nucleation of layers. Figure 7 (Plate LXIIT) 
shows (0001) monolayers spreading from a corner of a platelet, while fig. 8 
(Plate LXIII) illustrates nucleation from a contact with another crystal. 

Usually, the pyramidal and prism faces of all the hexagonal metals 
are coarsely striated, but some have high crystallographic quality. Even 
with phase contrast illumination, growth steps cannot be detected on 
these. This may mean either that they are absent or that they are too 
‘small to be seen with the present equipment. 
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§5. GrowTH PaTTERNS ON CADMIUM CRYSTALS 


These surfaces are very similar to those of magnesium crystals. Growth 
layers, probably about 5 A in thickness, whose visibility is slightly less 
than that of monolayers on magnesium, can be detected under careful 
phase contrast examination. Growth hills, centred on screw dislocations 
or nucleated from corners and edges of the platelets have been observed. 
Figure 9 (Plate LXIII) shows a spiral growth hill on a cadmium (0001) 
surface. 

§6. GrowrH PaTTeRNS oN Zinc CRYSTALS 
Zinc Platelets Grown Directly from Vapour 

The basal surfaces of the crystals grown in vacuo or in argon have 
similar quality except that the dark speckle which is characteristic of all 
zine growth surfaces is slightly denser when the crystals are grown 
in vacuo. The speckle, probably a coarse overgrowth of zine oxide, can 
be intensified by steaming the surface of a freshly grown crystal. It is 
often concentrated along the edges of growth layers which appear as 
dark irregular lines on a brighter background under phase contrast 
illumination. This ‘ decoration ’ of growth steps on zinc is well illustrated 
by fig. 10 (Plate LXIV). This shows monolayers which are probably 
only 5a thick; no measurements have been made. 

Although a large number of crystals have been examined and growth 
hills have been observed on many of these, there is no evidence for spiral 
growth of zinc. Usually, as illustrated in figs. 10, 11 (Plate LXIV), 
the step-lines form closed loops with a speck of material at the centre. 
Unlike magnesium and cadmium, there are several growth hills on each 
basal surface of zinc. This is a form of crystal growth from a low degree 
of supersaturation in which dislocations do not play a necessary role. 
The termination of growth steps in the surface indicates the presence of 
screw dislocations in the crystals but they are completely dominated by 
the layers from the two-dimensional growth centres. ; 

It is well-known that thin films of crystalline zinc oxide grow on zine 
metal with the basal planes of the two similar hexagonal structures in 
contact and with the [1010] directions coincident. There is, however, 
some doubt about the exact nature of the structure of the layers of zine 
oxide immediately adjacent to the zinc substrate. arly electron 
diffraction experiments (Finch and Quarrell 1933, 1934) indicated that 
these layers are deformed to account for the 20° misfit between the two 
lattices. This is the so-called pseudomorphic growth of zine oxide on 
zinc which has been discussed in connection with some theories of 
eryeualliie overgrowth (Frank and van der Merwe 1949, Mott and 
prey es Pots recent attempts (Raether 1950, Lucas 1951) 
i en eat Re he a techniques have failed to detect 
ie ane e structure of thin films of zine oxide, 
Rin ee one ick, growing on a basal plane of zine, appears to 

; a hexagonal network of zinc atoms in the oxide lattice 


with a= 3-22 A grows on a hexagonal network of zinc atoms in the metal 
lattice with a=2-67 4. 
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It seems likely that the growth layers observed on the basal surfaces 
are monolayers of zinc nucleated from small crystalline overgrowths of 
zinc oxide. The stock of subliming metal is a possible source of this 
oxide, since it is known that zinc oxide sublimes fairly easily if it is heated 
with zinc metal (Morse and White 1889). A zinc oxide crystal can serve 
as a nucleus for monolayers of zinc if it has stable size and shape and 
ensures coincidence of the hexagonal networks after a ‘ vernier repeat ’ 
(for example, a repeat of 35 zinc atoms in the metal substrate and 29 zine 
atoms in the oxide overgrowth along a common [1010] direction). 
As soon as the nuclei are exhausted, others begin to operate at other 
parts of the surface. It follows that a quantity of crystalline zinc oxide 
must be incorporated randomly throughout the metal crystal. Since 
the c-spacing of zinc oxide exceeds that of zinc by 6%, one zinc oxide 
crystal can initiate at most 17 successive layers of zinc without stacking 
faults: but if, as it appears, several zinc oxide crystals are active as 
growth catalysts at the same time, it is possible for correct growth to 
continue, merely producing periodic incoherence between the structures 
of the oxide crystals and the metal. Any place where an oxide crystal 
is buried in the lattice will, in general, be a centre of strain. 


Zine Crystals Grown by the Slow Cooling of Free Surfaces of Droplets 

The polycrystalline beads have a large number of randomly oriented 
facets. The appearance of hexagonal evaporation pits on most of these 
surfaces has identified them as basal (0001) habit faces of zinc. Some 
coarsely striated surfaces are obviously prism faces. 

Generally, the surfaces are similar in appearance to those of the zinc — 
platelets grown directly from the vapour. They have the characteristic 
speckle which often decorates the edges of monolayers. The phase contrast 
image of the steps shows that they form stepped valleys in the surface, 
probably due to evaporation. Figure 12 (Plate LXIV) shows that 
the step-lines form closed loops around dark hexagonal or triangular 
features ; these latter are probably polyhedral crystals of zine oxide. 
Again, zinc evaporates from the basal surface by a two-dimensional 
process, catalyzed by small crystals of zine oxide embodied in the metal. 
Another striking feature of the surfaces is the high density of small 
hexagonal or triangular pits with only one or two steps. The depths of 
these steps vary between about 5 and 1004. There are usually one or 
more large black clumps of material, also hexagonal or triangular in 
shape, at the centre of the pits (see fig. 13 (Plate LXIV) which is a 
phase contrast micrograph of a typical surface). It is thought that these 
are evaporation pits with large polyhedral crystals of zinc oxide at the 
centre. 


§7. GENERAL DiscUSSION OF THE GROWTH PATTERNS ON THE 
HexaconaL MBrALs 


Of the three metals examined, screw dislocations are responsible for the 
growth of magnesium and cadmium crystals from their vapours, but zine 
provides an interesting exception. The two-dimensional nucleation of 
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monolayers of zine is probably catalyzed by crystalline overgrowths of 
zinc oxide. The observations of growth markings on crystals grown 
directly from the vapour, and evaporation markings on the polycrystalline 
beads show that zinc oxide crystals must be incorporated during growth. 
It is interesting to note that the catalysis is frequently observed only 
on zine surfaces, and that zinc is the only hexagonal metal having an 
oxide with hexagonal structure ; magnesium and cadmium oxides are 
both cubic. It is difficult to understand why the growth of magnesium 
is not catalyzed by magnesium oxide more often, if the oxide grows 
with its (111) plane in contact with the (0001) plane of the magnesium 
substrate. The supply of oxide available during growth is probably very 
important. 

Usually, there are only one or two screw dislocations with Burgers 
vectors along [0001] directions in each crystal grown from the vapour. 
This is expected under the conditions of slow growth. Several crystals 
have crystallographically perfect basal surfaces but cannot be regarded 
as perfect crystals because they may contain edge dislocations or screw 
dislocations with no [0001] component. 


§8. Possrpiniry oF PLastic DEFORMATION EXPERIMENTS WITH 
THE HexacgonaL METALS 


Although some deformation markings which could represent slip lines 
with [112] directions in the basal surfaces of these metals have been 
observed, it is not easy to demonstrate that movement of screw 
dislocations could be responsible for these. So far, no movement of 
screw dislocations has been observed under very crudely applied stresses, 
but this is probably because other forms of deformation take place more 
easily. The most common of these is bending of the basal planes by a 
small angle, only a few minutes of are according to goniometric measure- 
ments. The ridge line in the basal surface, a diffuse line under phase 
contrast illumination, is sometimes approximately normal and sometimes 
approximately parallel to an edge of the hexagonal crystals; that is 
bending is along either [101] or [121] directions in the basal plane. 
Figures 1, 4, 5 illustrate these ridges ; a theory of phase contrast imaging 
(Hopkins, private communication) shows that the image of a sharp ridge 
should be a bright, diffuse band. Sometimes the bending is produced 
during growth, perhaps during the final cooling ; the step-lines are then 
observed to be kinked across the ridges. 
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EXPLANATION OF PLATES 


All photographs are taken under reflection, phase contrast illumination, x325. 
Fig. 1 
Showing a magnesium crystal containing a single screw dislocation. The 
step-line has been decorated by exposure to air and plasticine for 24 hours. 
Figs. 2 and 3 
Magnesium crystal with spiral growth hill on (0001) surface. 
Fig. 4 
Magnesium crystal with a growth hill centred on 13 similar screw dislocations. 
Fig. 5 
Magnesium crystal with a growth hill centred on 11 similar screw dislocations. 
Fig. 6 
Showing two-dimensional nucleation on a magnesium (0001) surface. Note 
the screw dislocations which are dominated by the growth layers on the 
right of the growth hill. 
Bin? 
Showing nucleation of monolayers from the corner of a magnesium crystal. 
Fig. 8 
Growth layers on magnesium nucleated from a contact with an adjacent crystal. 
Fig. 9 
A hexagonal platelet of cadmium with a single spiral growth hill. 
Fig. 10 
Showing a hexagonal platelet of zinc with several growth hills on the basal 
surface. The growth layers are thought to be initiated by small zine 
oxide crystals at the centres of the hills. 
Rig: 11 
Showing another zinc crystal with several growth hills with step-lines 
forming closed loops. 
Fig. 12 
Showing monolayers on a basal surface of zinc crystal grown by the slow cooling 
of a droplet. The monolayers are produced by evaporation which is 
initiated from the crystal of zinc oxide at the centre. 
Fig. 13 
An evaporation pit on a basal surface of zinc revealing a polyhedral 
crystal of oxide. 
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ABSTRACT 


The Doppler Effect has been observed for the gamma-rays of 4:5 and 
4-8 Mev emitted in cascade when alpha-particles of 958 kev are captured 
in 7Li. The gamma-rays are detected in a Nal(TI) crystal placed 
successively in forward and backward directions from the lithium 
target; the observed shift in energy between the two directions is 
(1-15-40-12)% and the theoretical first-order Doppler shift is 130%. 
It is remarked that these observations show the possibility of investigating 
certain gamma-ray life-times in the range 10~™ sec upwards. 


§1. INTRODUCTION 


Ir would be very surprising if gamma-rays did not show the Doppler 
shift. Such a shift has, indeed, been dramatically demonstrated for 
gamma-rays of 65 Mev from the decay of the 7° meson in flight (Bjorklund, 
Crandall, Moyer and York 1950). Doppler broadening of the 478 kev 
gamma-ray line from *Li* has been observed by Elliott and Bell (1949) 
and broadening and shifting of the same line has been observed by 
Rasmussen, Lauritsen and Lauritsen (1949). The success of the 
experiments of Moon (1951) on resonant scattering of gamma-rays 
imply the existence of a Doppler shift for the gamma-rays of 411 kev 
from !°8Hg. Du Mond, Lind and Watson (1949) and Hedgran and Lind 
(1951) have observed the Doppler broadening of the annihilation line 
from thermalized positrons. In none of these experiments, however, is 
the Doppler shift itself measured for an emitter of known velocity relative 
to the observer. The present note describes such measurements which 
demonstrate that the usual first-order Doppler formula is correct to 
within about 10% for gamma-rays of 4:8 and 4:5 Mev emitted from a 
body of speed 0-008 that of light. 


§2. GAMMA-RAY SOURCE . 


It was desired to put in evidence as large a Doppler shift as possible 
using a 1 Mv accelerator. Experiments such as those of Rasmussen et al. 
(1949) in which the gamma-ray is emitted from a nucleus recoiling from 
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the emission of a heavy particle suffer from the disadvantage that the 
velocity of the emitting body relative to the observer is compounded 
of that due to the capture of the bombarding particle and that due to 
subsequent emission of the product particle: unless coincidence 
techniques are used the velocity of the emitting body is not well defined. 
We therefore chose to use radiative capture in which case the momentum 
of the emitting body is that of the bombarding particle. The greatest 
shift may be obtained in the reaction ?H(py)®He, but this reaction is 
not resonant and so the stability of the bombarding proton energy is 
important. We chose to examine the gamma-rays from “Li(«y)!1B. In 
the range of alpha-particle energy accessible to us this reaction shows 
three resonances at 401, 819 and 958 kev (Bennett, Roys and Toppel 
1951). Of these the highest is strongest by a factor of about five. If, 
then, a thick target containing “Li is bombarded with alpha-particles of 
energy greater than 958 kev, the initial speed of the emitting body is 
well-defined as 0-00825c (if we neglect for the time being the effect of 
all but the highest resonance) and is in the direction of the bombarding 
beam of alpha-particles. 


§ 3. MEASUREMENTS 


The gamma-rays were examined with a 1-in. cube of Nal(Tl) and an 
E.M.I. multiplier type 6262. The resolution of this system was fair ; 
the pulse distribution from the 660 kev gamma-ray of !’Cs (strictly 
187Ba) had a full width at half maximum of 84% and a peak-to-valley 
ratio of rather better than 10: 1 for an uncollimated source. 

We found that the most frequent mode of de-excitation of the 9-27 Mev 
level in 11B formed at the 958 kev resonance is via the level at 4:46 Mev. 
We have worked with the combination of 4-81 and 4-46 Mev gamma-rays ; 
the 4-81 Mev gamma-ray is emitted first and the measured width of 6 kev 
of this resonance assures us that it is emitted in about 10-1* sec. It is 
probable that the cascade is of two electric dipoles and this in turn makes 
it probable that the life-time of the 4-46 Mev gamma-ray is about 107° 
(see Weisskopf 1951). Both these life-times are less than the slowing-down 
time of the 11B—of order 10~} sec. 

We found that the relative abundance of the two lines is a function of 
the angle of observation relative to the alpha-particle beam ; this demands 
that our angles of observation be symmetrically situated about the plane 
normal to the alpha-particle beam.* Our working angles were about 
48° and 132° relative to the alpha-particle beam. The centre of the crystal 
was 5-6cem from the target of LiOH and, when we take the average 
cos 6 over the crystal, the expected shift in energy between the two 
positions is 1-09%,. 

RS i ees ae See ek 2 ne 1s ee Se ee 

* All the states involved in this process appear to be sharp enough to 
prohibit the mixing of parities necessary to give a fore-aft asymmetry. 
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Owing to the smallness of the effect sought we could not trust the 
overall stability of the detecting system during the several counting 
periods of about # hour each needed to compile adequate statistics. 
Our procedure was to run for about 2 min on the gamma-rays from “1B 
then change the amplifier gain and run for a similar period and a similar 
number of counts on the gamma-rays from 17Cs, then back to the "B 
gamma-rays and so on. This gave a composite distribution of the two 
sets of gamma-rays; each distribution contained about 20 runs of each 
sort. The change in gain of the amplifier was arranged to bring the 
187Cy gamma-ray peak to a level comparable with that of the "'B peak 
but sufficiently far below it for the interesting portions of the two 
distributions to have a negligible effect on each other.* In this way 
the only element of the system on whose stability we had to rely to 
a high degree was the attenuator of the amplifier. In point of fact the 
stability of the whole system of crystal, multiplier, amplifier and kick- 
sorter changed by only 0-3% throughout the whole experiment of 6 hours. 
The kicksorter was of the type described by Hutchinson and Scarrott 
(1951), and we would like to express our thanks to Dr. G. W. Hutchinson 
for his kindness in loaning it to us. Six complete runs were performed, 
three at each angle, forming three pairs that were compared separately. 


§4. RESULTS 


Our procedure in analysing the results was to split each distribution 
into two parts : part A contained those channels that were predominantly 
filled with "Cs counts, and part B consisted of those channels above 
the 1°’Cs peak where the effect of the latter could be ignored. No use 
was made of the distribution between parts A and B. The distributions 
were plotted logarithmically both in numbers of counts per channel 
and in pulse height. Part A of a run at a forward angle was then 
compared with part A of its companion run at a backward angle; the 
two distributions were superposed on a light box and adjusted by eye 
until it was judged that the best fit had been attained. Each experimenter 
estimated that this fitting was made to about +0-10%f in pulse height 
and the independent judgments of the two experimenters agreed to 
within a comparable precision. Parts B of the same two distributions | 
were similarly and independently compared. In this way the fitting of 
parts A could not affect the fitting of parts B, and any possible changes 
in mean sensitivity throughout the two runs were eliminated. The 
fitting of parts B was to about -—.0-15%. This means that the overall 
estimated accuracy in each measurement of the Doppler shift (ratio of 
energies) is about +-0-18%. The three shifts obtained were 1:20, 1-22 
and 1-01%, giving a mean shift of (1:15--0-12)%. We may note that 
the spread between these independent determinations is of the order 
ne ee ee See 


- That part of the "'B distribution that lies beneath the 870s peak is rather flat. 
} Throughout this paper + means probable error. 
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of our estimated error in making the fit. As a check we compared 
parts A and B of distributions taken at the same angle and the 
independent fittings of the two parts agreed to within the expected 
limits. 

Figure 1 shows one of the three pairs of distributions; the scales 
have been adjusted according to the fitting in ordinate and abscissa 
of the ®’Cs peak, part A. It is seen that the 14B peaks (peak I is due 
to the 4-46 mev line, peak II to that at 4-81 Mev), part B, display a clear 
shift. The statistical errors may be estimated from the number of 
pulses involved.* Figure 2 shows an enlargement of the interesting 
region of part B. A further adjustment has been made between 
figs. 1 and 2 in ordinate only. 


S PER CHANNEL 


LOG COUNT 
Ke} 


es} 1-4 1-3 6 
LOG, PULSE HEIGHT IN VOLTS 


Pulse height distribution at 48° (filled circles) and 132° (open circles). Region 
A Ppprinepally due to 187Cs, region B to 141B—peak I due to 4:46 Mev 
line, peak II to 4-81 Mev line. 


§ 5. Discussion 


It is clear that the Doppler shift exists and that its magnitude accords 
roughly with expectation. A more detailed comparison of theory and 
experiment is complicated by several factors. One is that the peaks 
of part B do not represent the whole gamma-ray energy, but are due to 
pair creation in which both annihilation quanta escape from the crystal. 


ea . ; oa ee 
* There is little difference between the absolute ordinate scales for the two 
distributions. 
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In this case the expected shift would not be 1-0 % but rather 1-389. 
However, as is shown by a study of the 4-45 Mev line from *C, a strong 
subsidiary peak at higher energy corresponding to the capture in are 
crystal of one annihilation quantum exists, whose strength of about 50 Yo 
relative to the ‘ two-escape’ peak accords with expectation. For this 
peak alone the expected shift is 1:-22%, while for the relatively weak 
peak due to capture of both annihilation quanta the shift is 1-09%. 
When these complications are taken into account, as also the relatively 
small contributions from the lower resonances with smaller shifts, the 
expected shift for our method of fitting becomes 130%. This is rather 
larger, though not significantly so, than our experimental figure of 
(1:15+0-12)%, and we may consider that the Doppler effect for 
gamma-rays has been confirmed to within our experimental uncertainty 
of about 10%. 
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Enlargement of region B of fig. 1 with adjustment in ordinate only. 


It may be noted that the sense of the discrepancy is such as would 
be produced by a partial slowing down of the 11B before emitting the 
gamma-rays. However, as remarked above, we are certain that the 
4:81 Mev line at any rate is emitting in about 10-19 sec and that there 
are good reasons for supposing that the life-time of the 4-46 Mey line is 
only about 10-16 sec, while the slowing down time is of order 10-13 sec. 
Furthermore we see from the figures, particularly fig. 2 that the 4-46 mey 
line (peak I) is certainly shifted, and not noticeably less so than the 
4:81 Mev line (peak IT). 
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§6. APPLICATIONS OF THE METHOD 


In addition to verifying the Doppler shift this work shows the 
possibility of using the technique for the determination of life-times 
of certain excited nuclear states. In general we must suppose that 
coincidence techniques must be applied, together with provision for 
slowing down the recoiling nuclei in a determinate time. If this can be done 
it will become possible to measure certain gamma-ray life-times in the 
range 10~' sec upwards. We are at present engaged in applying this 
technique to the 6-14 Mev level of 180. 
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ABSTRACT 


Work on the energy expenditure per ion pair, V, by means of end- 
corrected ion chambers and proportional counters is extended to include 
ionization by Po «-particles. The gases argon, helium, hydrogen, nitro- 
gen, air, oxygen and methane are studied. Curves of ion current against 
percentage of argon admixed with the other gases are plotted and the 
ratio V gas/V argon determined in each case. Previous values of V 
obtained with electrons are discussed and small differences in the 
corresponding values of V in the electron and «-particle cases are noted. 


INTRODUCTION 


THE average energy expenditure per ion-pair, V, for electrons in various 
gases together with the ratio V gas/V argon have recently been investi- 
gated with proportional counters (Curran, Cockroft and Insch 1950, 
Valentine 1952). In the latter work a close agreement between the 
ratios V gas/V argon for slow electrons and for very energetic (340 Mev) 
protons (as obtained by Bakker and Segré 1951) was noted. In view 
of this agreement, it was decided to investigate the ionization of the 
same gases by Po «-particles using the same method as for electrons. 


EXPERIMENTAL RESULTS 


The method of using end-corrected cylinders as proportional counters 
and ionization chambers is immediately applicable to «-particles. The 
source was made by dipping a thick silver foil (area ~ em?) into a 
solution of Radium D+ E-+-F in hydrochloric acid. It was mounted on 
a probe running parallel to the central wire of the chamber and maintained 
at the voltage corresponding to its position in the field. The back of the 
foil was covered with ‘ Glyptol’ sealing paint so that all the «-particles 
were emitted towards the central wire. The energy distribution of the 
a-particles was examined by using the chamber as a proportional counter. 
The peak obtained was symmetrical except for a slight preponderance 
of low energy particles making the average energy of the pulse distribution 
97% of the peak energy. There are probably two causes for this energy 
loss, (i) finite thickness of the source, (ii) excitation of metastable states 
in argon; this energy would be lost in a proportional counter pulse, 


but detected in ion chamber operation. We deduce that the energy of 
eT eee ee i 
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the particles, so far as it affects ionization current, may be between 97 
4 100% of its known value and we adopt 98-5% as the most probable 
value. 

An absolute value of V for argon was found by determining the count 
rate of the source together with the saturation ion current produced by 
the «-particles in argon. After a careful examination of the various 
factors affecting the result the value V,—25-9 ev was obtained. Relative 
values of V were obtained from the ratios V gas/V,, which were found 
for the gases, helium, hydrogen, nitrogen, air, oxygen and methane, by 
measuring the ratio of the saturation ion current produced in the gas to 
that produced in argon. 

A further check on the relative values was obtained by plotting curves 
of ion current against percentage of argon admixed with the other gases. 
The numerical results are shown in table 1, columns 4 and 5, and the 
mixture curves in fig. 1. 


Table 1 


340 Mev 


Electrons Po «-particles 
protons 


Gas | Vinev | Vgas/Va | Vinev | Vgas/Va | Vgas/Va 


27-0 1-00 25-9 1-00 1-00 


32-5 ¢ 31-7 122 


38-0 37-0 1-43 


35'8 36-0 


DISCUSSION AND CONCLUSIONS | 


The absolute value of V,s—25-9 ev, obtained with Po «-particles, is 
in fair agreement with the value of 25-4 ev extensively used, but it is low 
compared with the more recent value of 28-3 ev obtained by Cranshaw 
and Harvey (1948). It should be noted, however, that their results on 
the variation of V, with «-particle energy are in marked constrast to 
those of Jesse, Forstat and Sadauskis (1950). A possible source of 
discrepancy is discussed by Hanna (1950). Jesse et al. measured the 
total ionization using a slow recorder while Cranshaw and Harvey 
measured only the immediately available ionization. The latter method. 
might miss delayed secondary electrons from metastable argon atoms and 
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so tend to give higher values for V. As the present work employed a 
total ionization method, we would expect the lower value. Indirect 
evidence in support of the value of V,4=25-9 ev is obtained by considera- 
tion of the values Vy,=36-0ev and V air=35-2 ev obtained using this 
value of V, together with the measured ratios. These derived values 
are in good agreement with the recent values of Vy,=—36-32-0-4 ev, 
V air—34-7-L0-5 ev obtained by Alder, Huber and Metzger (1947) and 
the value-of V air=35-1 ev, for Ra C’ «-particles, quoted by Gray (1944). 
The ratios V gas/V, are unlikely to be much in error and a value of 
V ,=28-3 ev would give Vy,=39-5 ev—a value higher than any obtained 
by recent workers. 


Fig. 1 
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Ratio of ion currents in gas mixtures. Argon with (a) methane, 
(6) helium, (c) hydrogen, (d) nitrogen. 


From table 1 it can be seen that the relative values obtained with 
a-particles are very close to those obtained with electrons except in the 
case of nitrogen and air. This may be due to an increase in Vy, for low 
energy electrons. From a consideration of previous results (Valentine 
1952), this increase is unlikely to take place for electrons of energy 
> 200 ev, in contradiction with the work of Pigge (1934) but in agreement 
with the results of Freund (1935) and Breunig (1936). 
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Values of V for the various gases were first obtained, using the above 
technique, by Curran eé al. (1950). Improvements in the method and 
the elimination of some sources of error (Valentine 1952) led to different 
values being observed in some cases. Table 1 summarizes the most 
recent values for electrons and «-particles obtained by the method and 
these values are considered to be more reliable than any of our earlier 
results. The sources of error have been carefully investigated and the 
values in table 1 are considered to be correct to about +2%. 


We should like to thank Professor P. I. Dee for his close interest in 
this work and his valued support. 
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ABSTRACT 


In this paper the Chapman—Enskog theory of non-uniform gases, as 
applied to the phenomenon of self-diffusion, has been utilized to evaluate 
the intermolecular force constants. The collision integrals occurring im 
the Chapman-Enskog formula for self-diffusion, as evaluated by 
Hirschfelder, Bird and Spotz by assuming the 12:6 power form of 
the Lennard-Jones potential function, have been used. The temperature 
variation of self-diffusion has been utilized to obtain the constant « of 
the Lennard-Jones function, which represents the maximum negative 
energy of interaction of two molecules. The intermolecular separation 
r) for zero interaction energy occurring in the same potential function 
has been calculated by actual substitution. The recent data on self- 
diffusion as reported by Winn, Winter and others have been used to 
obtain these constants and the values compared with those given by 
Hirschfelder et al. from viscosity and virial coefficient data. 


$1. INTRODUCTION 


In the kinetic theory of gases a quantity of great interest is the coefficient 
of self-diffusion, which is a measure of the rate at which a group of gas 
molecules diffuses into another group composed of identical molecules. 
Naturally it is not possible to follow the motion of individual molecules 
possessing no features distinguishing them from the other molecules of 
the gas, and hence no direct determination of the coefficient of self- 
diffusion could be made until recently. The only method available 
was an indirect one given by Lord Kelvin (1904) in which the coefficient 
of self-diffusion was evaluated by working with different triads of gases 
including the gas in question. Recently however, two important methods 
have been developed for measuring self-diffusion, in which, rigorously 
speaking, the inter-diffusion coefficient between two isotopes of the 
same substance is actually measured and this is taken to be equal to the 
self-diffusion coefficient of the gas. If the tracer isotope is a stable one, 
then the diffusion process is observed with a mass spectrometer ; while 
if it Is a radioactive isotope, the usual methods of detecting radioactivity 
provide a means of measuring the diffusion. In this way quite a large 


amount of data on self-diffusion at different temperatures have been 
ge 
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reported by a number of workers in recent years. It is the purpose of 
this paper to calculate the intermolecular force constants by utilizing 
these data and comparing the results thus obtained with previous 
determinations by other methods which are however for naturally 
occurring isotopic mixtures. 


§2. FoRMULAE AND METHOD oF CALCULATION 


The inter-diffusion of gases in general depends on the force fields of 
the diffusing molecules. Expressions for the coefficient of diffusion in 
terms of the interaction potential between the molecules have been 
derived by Chapman and Enskog and are as follows (Chapman and 
Cowling 1939) :— 

The first approximation to the coefficient of diffusion is 


(Os eo Ly AN ue ees wee ated avis (LY 
and the second approximation is 
[Dy2]e=[P19],/(1— 4), Soe oe ace See aD (2) 
where 
Pinot Pony te 
Ti Sag a trea 12 
( Q1%12+QeN1+ ie 
and 


P,=M E(u) 
P,,.=3(M,—M,)?+4M,M,A 
1= ME) (u)1 (6M? + 5M ,?—4M °B+ 8M M.A) 
Q12=3(M,—M,)*(5—4B)+4M,M,A(11 —4B)+ 2H? MM 4/(H4)1(H2)1 
with similar relations for P, and Q,; and 


My=M1+ Me, M,=m,/mM, M,=m,/mM, 
Np=N1 +N; Nyy=N4/N, Nog=No/Nq; 
Ny2=N4/Ng, No1—=No/N4, 


where the m’s are the masses of the gas molecules and n’s their number 
densities. (,), and (yw), are the first approximations to the viscosities 
of the two components, given by 5k7/82,° in each case ; the quantities 
A, B, C and E£ are functions of the collision integrals given by 

A= 2,752, ; B=[52,— 2,M52,™, 

C6=29,P/52,%, and H#=kT/8M,M,2,, 
and the various collision integrals 2,” refer to the interaction of a 
molecule of component 1 with one of component 2. 4 is a small 
correction so that [D,,.]. rarely differs from [D,,], by more than 1%, 
which is within the limits of experimental error. We shall therefore 
take [D,,], to be equal to the experimentally determined values of the 
diffusion coefficient. The first approximation [D,,], depends only on 
encounters between molecules of opposite type and thus to a first 
approximation, encounters between like molecules do not affect diffusion. 
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The collision integrals Q,(” occurring in the above formulae depend 
upon the form of the intermolecular potential function. As already 
discussed by Buckingham (1938) and others, the interaction energy 
between two molecules for a separation 7 is best represented by the 
Lennard-Jones equation of the form 

E(r)=4e[(ro/r)2#—(ro/r)1, - . - - + + + (83) 
where 7) is the separation for which energy of interaction is zero and 
« is the energy difference between the separated molecules and the 
molecules in the configuration for which they have the maximum negative 
energy. Hirschfelder, Bird and Spotz (1948) have carried out the 
numerical evaluation of the collision integrals 2, for this form of the 
interaction potential function and have tabulated the values of the 
functions W'(n; x) which are related to the 2, of Chapman and 
Cowling by 

QM =12(2akTm,/mym,)12W'(n ; 2), 
where «=e/kT’. 
Hirschfelder et al. on the basis of the Chapman-Enskog theory give 
for the diffusion coefficient the expression 
0-000929167[(M,+M,)/M,M,]¥? 
OE ae ee (4) 
Pr 2) Wi; kT /ey2) 
in which D,, is the coefficient of diffusion in cm? sec~1, p is the pressure 
in atmospheres and M,, M, represent, here, the molecular weights, 
and W1(1; kT'/e,.) is a slowly varying function of kT'/<e,.. Equation (4) 
is readily deduced from eqn. (1). It is obvious from the above equation 
that 7?/?/[D,,.], will obey the same variation with temperature as W,1 
provided 1. is assumed to be constant. It is therefore possible to evaluate 
€,2 (provided it is assumed to be constant) from the observed temperature 
variation of 7?’?/[Dyolexp oY T?[D,.], and from this value of ¢,., the 
value of 7/19) the zero-energy intermolecular separation for molecules 1 
and 2 can be found by actual substitution in eqn. (4). 

This procedure is adopted below for oxygen, nitrogen, argon, neon, 
carbon dioxide and methane, for which the data for variation of self- 
diffusion coefficient with temperature are given by Winn (1947, 1948, 
and 1950), Winter (1951), Hutchinson (1949), Groth and Sussner (1944) 
and Harteck and Schmidt (1933). The data of Winn give D,, which is 
calculated from the relation (as shown by Hutchinson) 


D1,=[2M,/(M,+M,)}""[Dyslexp; 

where M, and M, are the molecular weights of the two isotopes. His 
data have therefore been reconverted to give [Djs]exp for the diffusion 
of one isotope into the other which is the quantity actually measured. 
by him. In order to determine ¢,,, the following procedure was adopted : 

Hirschfelder’s tables were utilized to give a plot of W,} against kT'/e, 
from which a second graph was plotted for [W,+],/[W,1], against various 
initial values of [W,1], for a given ratio of [AT'/e]o/[kT'/e],, the suffixes 1, 2 
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outside the square brackets refer to temperatures 7, and 7',. Such 
graphs were plotted for the values 1-5, 1-8, 2-0, 2-5 and 3 of the ratio 
[kT'/e],/[kT/e], and are given in fig. 1. Then a graph of T?/[Dyslexp WAS 
plotted against 7’ (fig. 2) which served to check the accuracy of the 
experimental data and further gave us a number of additional points 


Fig. 1 


[wil,/fwi] 
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for calculation. From the 7?/?/[D,s]exp versus 7’ graph, values for 
[7'/?/(D1)explo/[T?/?/(D12)exp]i for the various ratios of kT Je were calculated. 
at a number of initial temperatures from the graph and for each calculated 
value of [73!?/(D42)exp]o/ [T??/(D12)expl, the corresponding value of W,* 
was read from the appropriate ratio graph of [W,1],/[W,1],. From this 
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value of W,} the corresponding initial value of kT',/« was read from W,+* 
versus kT/e graph. Knowing the initial temperature which we took for 
[Dyolexp 01 T?#/(Diz)exp graph the constant « was obtained by simple 
substitution. In this way a number of values were obtained for €), and 
the arithmetic mean calculated. 

Fig. 2 
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§3. RESULTS 


The observed data, mostly of Winn (1950), are plotted in fig. 2 for 
argon, neon, oxygen, nitrogen, carbon dioxide and methane. For argon 
the ratio [T?/?/(D4)explo/[T'*?/(D12)expl1 calculated from 7',/7;=2,/%,=1-5 
at initial temperatures 192° K and 204° K comes out to be 0-869 and 0-875 
respectively ; and for x,/2,=2 at initial temperature 90° K, it comes 
out to be 0-711; for w,/a,=2-5 at initial temperatures 90° K and 105° kK 
it comes out to be 0-646 and 0-667. Similarly for x,/7,—3 at temperatures 
90° kK and 105° K it comes out to be 0-611 and 0-633. These being read 
a eh respective [W,1],/[W,"], graphs yield for the ratio 1-5 the values 
ob W,* as 0-580 and 0-562; for the ratio 2 the value is 0-890 ; for the ratio 
2-5 the values for W,1 are 0-882 and 0-80 and for the ratio 3 the values 
are 0-865 and 0-785. These by reference to the W,! versus kT/e graph 
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gave the values for ¢,,/k as 117-1, 118-5, 130-4, 131-3, 128-6, 127-3 and 
125-0 giving the mean value 125-5. Using these values of ¢,, the 
corresponding 79:15) values as calculated from eqn. (4) came out to be 3-435, 


Table 1 


Argon Neon 


Temp. Range  «/k roin A Temp. Range e/k 
re) K ° K 


90-180 130-4 3°331 

90-225 128-6 3°346 75-187-5 

90-270 125-0 3°378 75-225 
105-262:5 1313 3371 84-168 


105-315 127°3 3-403 90-180 
192-288 118°5 3°423 90-270 
204-306 UM eal 3-435 


Mean 125-5 3°384 Mean 


Oxygen . Nitrogen 
Temp. Range e/k Temp. Range e/k 
Ue or 


114-342 129-5 75—-187°5 101-4 
135-337°5 110-6 75-225 100-0 
150-270 130-4 90-225 112-6 


150-300 118-1 90-270 103°5 
165-297 117-9 120-300 102-6 
165-330 108-6 150-300 94:3 


Mean 119-2 Mean 102-4 


Carbon dioxide Methane 


Temp. Range e/k ro in A Temp. Range elk 
fe) ike fo] K 


180-360 145-2 4-135 196-352-8 176-6 
195-351 147-7 4-124 220-330 157-2 
210-315 156-7 4-063 224-336 138-2 
225-337°5 146-1 4-122 228-342 144-5 


Mean 148-9 4-111 Mean 154-1 
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3-423, 3-331, 3-371, 3-346, 3-403 and 3-378, giving the mean value as 
3.384 in A.u. A similar procedure was adopted for other gases also. 

The €,5 and roy) values thus obtained are given in table I for the gases 
investigated, together with the temperature ranges for which they have 
been calculated. These results are conipared in table 2 with those of 
Hirschfelder and others from data on viscosity and virial coefficients. 


Table 2 


From From From From Viscosity} 


_ self-diff. Viscosity* | virial coeff.* 


Gases Low temp.| High temp. 


lg e/k ro e/k to to e/k to 


Argon 3-384 |124-0 | 3-418 ‘0 | 3-483 | 126-2 | 3-414 
Neon -2 | 2-659 | 35:7 | 2-80 35-7 | 2°74 ‘0 | 2-769 | — — 
Oxygen 3-260 1113-2 | 3-433 ‘7 | 3-435 | 134°9 | 3-325 
Nitrogen 3-478 | 91:46] 3-681 | 95-9 ‘8 | 3-687 | 117-1 | 3-537 
CoO, 4-111 |190-0 | 3-996 
Methane 3°648 |1386:5 | 3-822 


* Hirschfelder e¢ al. (1948). + Srivastava and Madan (1952). 


§4. Discussion oF RESULTS 


For an exact evaluation of « the accuracy required in the measurement 
of D,, is too high to be realized experimentally. The experimental error 
in D,, is usually of the order +1% which would produce an error of 
about +15% in the value of «. Thus the value of « as determined from 
measurements of self-diffusion can be taken to be correct to +15%. 

For neon, it will be seen from table 2 that the value of « from self- 
diffusion is higher than that calculated from viscosity. It may be 
remarked that the neon molecule is known to be much harder and it is 
possible that this discrepancy might be due to that cause. 

For nitrogen and carbon dioxide the observed deviations may be due 
to insufficient data at low temperatures. 

It is seen from table 1 that in the small range of temperatures over 
which diffusion measurements have been made there is no systematic 
change in the values of ¢,. and 794.) as the temperature is progressively 
increased. This shows that an equation of the form (4), with the 
intermolecular potential as given by (3), is fully confirmed by the 
existing data and is a simultaneous verification of the Chapman—Enskog 
theory and the Lennard-Jones form of the intermolecular potential 
function. 

It may be remarked that values of ¢« from viscosity and virial 
coefficients refer to the normal gases and may thus be taken to apply to 
the most abundant isotope (or in the case of neon to the normal mixture). 
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Theoretical considerations seem to show that it is unlikely that there will 
be any measurable difference between the coefficient of self-diffusion 
for the most abundant isotope, and the coefficient of inter-diffusion for 
that isotope and any other. It is therefore not possible to say whether 
the slight differences found here between the values calculated from 
self-diffusion data and from viscosity data are due to experimental errors 
or due to an actual difference between the two isotopes. 

The foregoing calculations are based on the assumption that both « 
and 7) are constants. The consequences of their variation with 
temperature are being investigated and will be reported later. 


We record our thanks to the Uttar Pradesh Government for the 
award of a research grant. 
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ABSTRACT 


Surface measurements are made upon kink bands, twin bands and 
bend planes produced in specimens cleaved from single crystals. The 
microscopic structure of kink bands is illustrated, while the method used 
to determine the twinning elements indicates relations in the shape 
of the twinned region. Polygonization structures produced by 
inhomogeneous deformation followed by annealing are illustrated inter- 
ferometrically and related to growth features found in metal crystals. 
In the case of hexagonal close-packed crystals, deformed by glide, a 
type of deformation band structure of the active glide plane is illustrated. _ 


$1. INTRODUCTION 


THE methods of interferometry (Tolansky 1946), can be applied to the 
measurement of the topography of deformation structures found upon 
metal surfaces in those cases where a high degree of perfection of the 
surface can be attained initially. In the case of the metals which possess 
good cleavage properties the surface revealed by cleaving at liquid air 
temperatures fulfills this requirement well; furthermore this surface is 
of particular interest in plastic deformation studies since the plane of 
good cleavage is also, in most cases, the prominent glide plane. A recent 
review of the deformation processes of zinc crystals (Jillson 1950), has 
drawn attention to the variety of markings found upon the cleavage 
surfaces of that particular metal. An introductory survey of the structure 
of markings found on the cleavage planes of Zn, Mg, Bi and Sb crystals 
is made in the present paper, the quantitative methods of interferometry 
being of particular use in classifying the markings observed on the 
cleavage surface after the crystals had been subjected to plastic 
deformation experiments. 


§2. GrowtH FEATURES 


Single crystals of the cleavage metals were grown from relatively 
high purity material (Zn 99-997%, Mg, Bi and Sb 99-99%), by the usual 
travelling furnace method with the modification that the growing crystal 
was free from undue restraints from its container. When growing crystals 


of metals which twin readily, this is essential in order to avoid large 
Ee eee 
* Communicated by the Author. 
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regions of the resultant crystal having a twin orientation. The well known 
imperfection of metal single crystals grown from the melt, that the 
cleavage surface is frequently not perfectly plane but displays boundaries 
between areas of slightly differing inclination, is suited to interferometric 
illustration. Figure 1 (Plate LXV) shows interference fringes crossing such 
an ‘ optical mosaic ’ boundary in the case of a Zn cleavage. The inclina- 
tion between the two basal areas is 1° 2’, a typical value for crystals which 
display only one or two mosaic boundaries on their cleavage face. 
Survey of the entire cleavage shows that the boundaries although 
smoothly curved tend over most of their length to lie at right angles to 
a possible glide direction. The area illustrated in fig. 1 is typical in this 
respect, as can be seen from the position of the boundary with respect 
to the twin band and associated bend plane produced in the crystal 
mechanically. Furthermore, in specimens displaying a number of 
boundaries the consequent tilting of the various cleavage areas forms a 
cleavage which is either completely concave or completely convex. 

The growth conditions which affect the presence of optical mosaic 
boundaries have been extensively discussed (Holden, A. N. 1950), but the 
feature is of interest to the plastic deformation experiments to be 
described, in that the latter can lead to the production of such boundaries 
within the crystal in a regular manner. 


§3. Twin Banps AND BEND PLANES 


The cleavage of plate-shaped specimens from a single crystal rod 
results in cleavage surfaces which display the traces of twinning move- 
ments at various stages of development and which owe their origin to 
the local cleavage impact. Many of the twinning movements do not 
progress very far into the crystal and display on the cleavage surface 
the typical lenticular trace. The accommodation of these limited twin 
regions within the original crystal is brought about by the sharp bending 
of the basal planes of the crystal about the position of termination of 
the twinning movement on the twinning plane K,. Figure 3 
shows diagrammatically the shape of such a twinned region as 
determined by successive sectioning of a specimen while fig. 2 (Plate LXV) 
illustrates interference fringes contouring a section of the twin trace 
on the cleavage plane in the region XY. It can be seen by following 
a fringe of given order through the section that the amount of bending 
at the bend plane is that which is necessary to accommodate the twinned 
region within the larger restraining matrix. Under less restrained 
conditions of twinning, when for example a plate specimen is subject to 
an external shear in an appropriate direction, bend planes are not seen 
since the twinned regions develop to such a degree as to pass completely 
through the specimen and occupy a considerable volume. Such twinned 
regions are illustrated in fig. 4 (Plate LXV), the surface being the (II) 
cleavage of Bi. It will be observed that the interference fringes give a 
precise measurement of the angles produced on the crystal surface by the 
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twinning movement and in the case of the crystallographic cleavage 
surfaces it is only necessary to measure the angle made by the twin band 
on the cleavage surface to determine the shear, s, involved in the twinning 
movement. For non-cleavage metals the method of Mathewson (1928) 
for determining the twinning plane K, can be extended to include the 
determination of the direction and amount of shear in twinning by the 
use of interferometric goniometry such as is illustrated in fig. 4 (Plate LXV) 
(Holden, J. 1951). 

From fig. 2 (Plate LXV), using the notation of fig. 3 it is seen that 
61 =bB where 9 and b are the angles produced on the cleavage surface 
by the twin band and bend plane respectively. T' is related to T,, the 


Fig. 3 


Cleavage 
Plane 


Shape of twinned region. 


development of the twin in a direction normal to the twin plane K, 
by the relation 7,,=T7 sin ¢ and B is related to the development of ihe 
twin in the direction of shear y, by the equation B=D cos @; being 
the angle between the plane K, and the cleavage plane. The extent of 
the twin region in the direction normal to n and 7, is seen as the trace 
of the twin in the cleavage plane, this being the direction of greatest 
development, the twinning movement spreading at a very early stage 
of its development to the limits of the crystal in this direction, It F 
thus seen that in the early stages of the twin development the ratio of 
the extent of the twin region in the directions 9, and mn is given b 

(D/T,,)=6/(6 cos¢ sind). In the case of twinning in Zn, @=4° oo 


Re 
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p=47° 2’, while values of 6 can be directly measured for a range of 
restricted twinning movements such as are found during the cleavage 
of plate specimens. The ratio is found to be constant and characteristic 
of the particular metal crystal for twinning movements produced under 
similar conditions. In the case of the cleavage of Zn plates at liquid 
air temperatures values of b=44+-2' are found giving (D/T,,)=11-0-5, 
while for Mg under similar conditions (D/T',)=13-.0-5. The sharply 
positioned bend plane would suggest that the development of a twinning 
movement within a crystal in the directions n and y, is dependent upon 
a limitation of the maximum angle of bending which the bend planes 
can attain in order to accommodate the twinned region. This is of interest 
in connection with the properties of individual boundaries as produced 
in the polygonization of hexagonal crystals. As will be described in a 
following section such polygon boundaries are subject to an upper limit 
of approximately 1° to the change of orientation which can take place 
at any one boundary. Figure 5 (Plate LXV) illustrates the termination 
of a twin region as indicated by its trace on the cleavage plane; the 
narrow twin band is out of focus in order to achieve adequate contrast 
of the bend plane running alongside and terminating with the twin 
band. The close relation between the bend plane and the twin band 
traces can be observed when a local indent of the cleavage surface causes 
the collapse of part of the twin region or where twin regions react with 
each other. In these cases the position of the bend plane indicates that 
the proportions of the twinned region in the directions n and 7, fall within 
the limits previously noted. 


§4. Kink Banps 


Cleavage of a Zn crystal rod by the gentle impact of a sharp edge 
depends upon orienting the cleavage blade parallel to the basal plane. 
Thus the crystal is in a position in which it may suffer compressive 
stress along the basal planes. This condition has been shown on a 
macroscopic scale to lead to the formation of kink bands (Orowan 1942). 
The typical bands formed on a microscopic scale during the cleavage of 
a Zn crystal at liquid air temperatures are shown in fig. 6 (Plate LXVI). 
The trace of these bands is at right angles to a glide direction and it is 
seen that the individual markings collect together into relatively broad 
bands. This is in contrast to the fine scale markings seen on the cleavage 
faces of crystals which only twin under the compressive stresses associated 
with cleavage. In the case of Bi, for example, fine markings on the same 
scale as kink bands are observed, however, they are more evenly spaced 
and also run parallel to the trace of twin planes. The characteristic 
topography of a collection of kink bands in Zn is illustrated in fig. 7 
(Plate LXVI) in the case of a specimen in which by choice of appropriate 
cleavage direction only one set of bands have been produced. It is seen 
that the broad bands are places where individual kink movements of 
opposite arrangement have piled up. The symmetry of the piling is 
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marked in fig. 7 (Plate LX VI) but it is always found for all broad bands that 
the location of the band is about the meeting of individual kink bands of 
opposite arrangement. Figure 8 (Plate LXVI) illustrates a particularly 
simple case of the resultant position of individual kink bands ; the local kink 
at the crystal surface is most precise and the structure performs under 
suitable illumination conditions as a diffraction grating of the appropriate 
number of ‘slits’ and spacing. ‘The change of height of the crystal 
surface at an individual kink band can be measured interferometrically, 
but the measurement of the distance along the surface over which this. 
height change takes place calls for careful phase microscopy. The 
dimensions of the individual kinks of fig. 8 (Plate LX VI) are found to be as 
follows: change of level of the surface at the kink ; 9001004; width 
of the kink 3-8--0-4; distance between pairs of opposite kinks 45-+-5 p. 
The kink structure was produced on a cleavage plate approximately 
1 mm in thickness at liquid air temperatures. The surface angle for the 
isolated kink bands is thus 1-4--0-3° and this order of magnitude is again 
found for the individual kinks within the broad piled up bands. 

The cleavage surface of Zn crystals shows kinking movements freely 
when the crystals are suitably compressed (Jillson 1950), or indented. 
Individual kink movements can also be studied microscopically on 
Zn cleavages since the termination of a restricted twinning movement 
is usually associated with the generation of groups of individual kink 
bands. Precise measurements upon the dimensions of these kink bands 
in their early stages of development are at present being attempted. 


§5. INHOMOGENEOUS DEFORMATION 

When a plate cleavage specimen free from surface markings (in 
particular the tilt boundaries of growth) is bent to a radius of a few 
centimetres about an axis lying in the cleavage plane and at right angles ~ 
to a glide direction, then annealed for, say, 12 hours at 100—-150° co below 
the melting point of the metal, the originally curved cleavage surface 
becomes faceted, the boundaries of the facets running accurately at 
right angles to the appropriate glide direction. Figure 9 (Plate LX VII) 
illustrates the facets on a Zn crystal plate treated in this manner, the 
radius of bending being 3-2 cm and the annealing time 10 hours at 300° c. 
The interference fringes indicate that the facets are inclined to each 
other at angles in the range 34’ to 42’. Identical results are obtained if 
an entire single crystal rod is bent and annealed and then the cleavage 
specimens obtained. The cleavage thus faithfully follows the slight 
inclinations of the cleavage plane. By cleaving successive thin plates 
off the crystal rod the three dimensional picture of the facet system can 
be built up. In the case of the bending axis described the facet system 
changes only gradually over the length of the crystal rod, some facets 
prominent in one section becoming less prominent and being replaced 
by others as the sectioning moves over lengths of some centimetres. 
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Should the axis of bending while still being in the cleavage plane not 
coincide with a direction at right angles to a glide direction then the 
facets at any cleavage appear bounded by polygons, as illustrated by 
the micrograph of fig. 10 (Plate LXVII). As the bending of the specimens 
becomes more severe it is found that the number of facets increases, since 
there appears to be a limit to the angle possible between the facets 
produced by the annealing, which for the relatively pure Zn crystals is 
about 1°. This upper limit can be demonstrated in the case of the 
local curvatures produced by allowing a steel ball bearing to fall upon 
the cleavage plate. On the opposite face of the plate to that receiving 
the impact is raised a shallow dome ; by choice of suitable balls a wide 
range of local curvature can be produced and accurately measured 
interferometrically. When such a specimen is annealed on a hot stage 
microscope the first effect observed is the appearance of a number of 
concentric rings at the position of the dome ; if the annealing temperature 
is relatively high the rings soon become quite separate and can be 
counted. The rings mark areas of slightly changing tilt which can be 
directly measured at any particular ring or the number of rings can be 
counted and in conjunction with the known initial dome curvature 
an average value of the change of inclination at each bounding ring 
arrived at. For Zn, measurements made upon a wide range of local 
curvatures indicate that the change of inclination at a boundary does 
not exceed 1° 5’. This value refers to material of 99-997°% purity and 
annealing temperatures up to 350° c. Continued annealing of a specimen 
causes the rings to develop into hexagons, the sides appearing in positions 
such that their traces are at right angles to directions of glide. At the 
same time the boundaries move away from the position of the dome thus 
eliminating it as a feature. Progress of the hexagon boundaries through 
the crystal continues over a period of hours the distance between 
neighbouring boundaries increasing, eventually it is seen that highly 
symmetrical patterns of polygons are formed away from the original 
dome position in the directions of glide on the crystal face. In these 
directions the curvature being relieved is greatest about one glide 
direction and smaller and equal for the other two directions. In this 
movement of tilt boundaries any twin bend planes which may be produced 
in the crystal by external means enter into the composite patterns. By 
their method of production and properties the boundaries described are 
polygon boundaries (Lacombe 1948, Cahn 1948). The temperature of 
annealing could be decreased and the time of annealing increased within 
wide limits to produce similar structures. The traces of kink bands were 
apparently unaffected by such annealing processes. 

The manner in which the polygonization structures orient themselves 
suggests the origin of the optical mosaic phenomenon mentioned in 
connection with the growth of single crystals. It is observed in metal 
erystals grown under conditions favourable to the production of numerous 
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optical mosaic boundaries that cleavage of the crystals shows that the 
boundaries form on the cleavage surface a system of very irregular 
polygons. However, one system of sides tends to predominate and to 
lie in a direction at right angles to a glide direction. The limit to the 
maximum angle of inclination at a polygon boundary, as determined in 
the deformation experiments, also indicates an interpretation of optical 
mosaic structures in terms of polygonization, since in the case of good 
single crystals (crystals containing only one or two mosaic boundaries), 
the residual misorientation between the parts of the crystal after the 
growth process is generally close to the upper limit characteristic of 
polygonization (cf. fig. 1 (Plate LXV)). 


$6. GLIDE DEFORMATION 


Cleavage of a single crystal rod of Zn or Mg can be carried out after 
the crystal has been subject to tensile stress and caused to glide. 
A specimen some 20cm long being used in order to allow the cleavage 
to be made well away from the gripped ends. When the specimen is 
extended slowly at a moderately high temperature the cleavage face 
shows to the eye the appearance of widely spaced straight bands running 
completely across the face accurately at right angles to the active glide 
direction. Interference fringes contouring these bands show that each 
local band is due to the juxtaposition of equal and opposite local 
curvatures of the glide plane. Figure 11 (Plate LXV) shows the topography 
of one band in the case of a Zn crystal subject to 12% glide strain in 
a time of 30 minutes, the crystal being cleavaged at a random position 
well away from the gripped ends. The band illustrated was one of six 
evenly spaced on a cleavage face measuring approximately 18 mm in 
the active glide direction. The bands extended into the crystal and were 
normal to the glide direction as could be observed by cleaving the crystal 
rod into thin plates. This sectioning showed that the individual bands 
had their maximum curvatures at a given section and slowly faded 
away on each side being replaced by other bands coming to their 
maximum curvature. With increasing glide strain the number of bands 
was not found to greatly increase in the case of these hexagonal crystals, 
however, the curvature within the existing bands was found to build up 
in a regular manner. 

Deformation at room temperature at similar rates to the higher 
temperature deformations resulted in band formation of the same type 
but the bands were, in general, not so completely discrete and symmetrical 
and did not completely cross the crystal face. Unlike the higher 
temperature bands, which for a given glide strain were few in number 
and almost identical in curvature, the more rapid room temperature 
extensions to the same glide strain gave rise to more numerous smaller 
bands (that is in terms of their superficial extension on the glide plane in 
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the direction of glide), of a more diverse range of curvature. The hardening 
of a hexagonal single crystal rod can be rapidly attained by bending 
and unbending, for simplicity about an axis in the cleavage plane at 
right angles to a glide direction. In this case the cleavage face is found 
to be traversed by a close packed system of bands at right angles to the 
glide direction. The corrugated nature of the glide plane under these 
conditions is illustrated in fig. 12 (Plate LXVI). 

The large scale S-bands and the finer scale corrugations can be 
polygonized by annealing, since they consist essentially of local 
curvatures. Figure 13 (Plate LXVII) shows the detail of a polygonized 
S-band in a Zn crystal, the band having been made to be of particularly 
marked curvature (by subjecting the crystal to considerable glide strain), 
in order to show up in the micrograph. The maximum angle of tilt at 
a boundary formed by the polygonization of an S-band was found to be | 
within the limit already stated for polygonization structures of simpler 


type. 
$7. SUMMARY 


The consistent structures found upon the glide plane of metal crystals. 
subject to a range of plastic deformation experiments, and revealed by 
cleavage, substantiate the method as a means of preparing a metal 
surface for precise measurement. The methods of interferometry permit. 
illustration of the three dimensional structure of twinned regions, kink 
bands, polygon boundaries and glide deformation bands and also 
lead to measurements bearing upon theoretical concepts of crystal 
imperfections. Further measurements upon fine scale structures found 
upon the cleavage surfaces of the metals used in the present survey 
and upon other metals are in course of publication. 
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EXPLANATION OF PLATES. 


Puate LXV. 

Fig. 1 
Optical mosaic boundary, Zn. (x 120) 

Fig. 2 
Twin band and bend plane, Zn. (x 100) 
Fig. 4 

Large twinned regions, Bi. (x 40) 
Fig. 5 

Termination of twinned region, Mg. (x40) 


Fig. 11 


Topography of glide deformation band; the vertical represents the glide 
direction, Zn. (X30) 


Puate LXVI. 
Figs. 6, 7, 8 
Kink bands, Zn. (x 220) 
Fig. 12 
Deformation bands, Zn. (X30) 


Puate LXVII. 


Fig. 9 
Polygonization boundaries ; the horizontal represents the glide direction, Zn. 
(x50) 
Fig. 10 
Polygonization boundaries, Zn. (30) 
Fig. 13 
Polygonization of glide deformation band, Zn. (x40) 
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ABSTRACT 


The momentum distributions in the centre of mass system for the charged 
secondary from the two reactions, V—-p*+7°+yv and V+—-p+-+vr-+4, 
have been calculated. The second of the above processes has been 
used to calculate the transverse momentum spectrum of the charged 
secondaries. This has been compared with the measured spectrum from 
charged V-particles and found to be consistent with it. It is not possible 
to distinguish between the two decay schemes by a study of the momentum 
distributions, but alternative methods of distinguishing between them 
are discussed. 


§1. INTRODUCTION 


In a detailed study of charged V-particles, Armenteros et al. (1952) 
have given evidence that the transverse momentum distribution of the 
decay products did not fit a single two body decay scheme. It is the 
‘purpose of this paper to calculate the transverse momentum distribution 
assuming a three body decay scheme and to compare this distribution 
with the V-particle data. Two possible decay schemes have been 
suggested, viz. 

Vt—+pt+qr+y, Te Ge a he See. 2 (1) 


Vtsputtv+y. 27 phigh iene a2) 


‘These two modes of decay will be considered in the following calculations. 


§2. CALCULATION OF THE TRANSVERSE MOMENTUM DISTRIBUTION 


The charged V-particles are observed to decay in flight and the trans- 
verse momentum distribution of the charged secondaries can be readily 
obtained from the experimental data. Because this distribution is 
invariant under a Lorentz transformation, it will be calculated in the 
centre of mass system of the decaying V*-particle and then compared 
with the distribution observed in the laboratory system. 
at a eM 2 ae tS 


* Communicated by Professor P. M. 8. Blackett, F.R.S. 
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If a three body disintegration scheme is assumed, the transverse 
momentum distribution of the decay products of charged V-particles 
can be analysed in terms of two factors. The first is that in the centre 
of mass system for the decay, the charged secondary will be emitted in 
any direction with equal probability. This angular distribution will 
cause a transverse momentum distribution for any given momentum of 
the charged secondary particle in the centre of mass system. The second 
factor is that, in a three body scheme, the momentum of the secondary 
particle will not be a constant, but will vary depending upon the mode of 
decay. 

As a result, the transverse momentum distribution must be expressed 
as the product of these two independent distributions. Because one 
does not know the momentum of the secondary in the centre of mass 
system from the experimental data, it is necessary to average the angular 
distribution over the momentum distribution in the centre of mass system 
to obtain a transverse momentum distribution which is directly com- 
parable with the experimental data. 


(a) The Angular Distribution 

Armenteros et al. (1951) have shown that for a given value of the 
momentum, p, in the centre of mass system, the transverse momentum, 
pr, is distributed according to the expression 


prdpr 
PV (P°—Pr")’ 
as a result of the random orientation of the direction of the charged 
secondary. This then is the angular distribution of pp, which must be 
averaged over the distribution of values taken by p. 


(3) 


(b) The Momentum Distribution in the Centre of Mass System 


Michel (1950) has given a general treatment of the interaction of four 
Fermions. In a later paper (1952) he has calculated the energy distri- 
bution for the y-meson in the decay scheme (2) above. We can obtain 
the differential momentum distribution from the generalized form of 
the energy distribution given in Michel’s first paper. As a basic simpli- 
fication to this expression, however, let us assume that the matrix element 
of the interaction is independent of the energies of the particles produced 
by the V*-particle disintegration. If this approximation is made, it 
is no longer possible to distinguish between Fermions and Bosons, and 
the form of the energy distribution is determined by the density of states 
available to the decay products in phase space. This approximation 
seems justified in that the general form of the distribution is preserved 
and the experimental data available for comparison is very meagre. ; 

Let us call g the coupling constant of the reaction, M the rest mass of 
the decaying particle, » the rest mass of the charged decay product, p 
the momentum of this charged secondary particle, and X the rest mass of 
one of the neutral secondary particles. The other neutral decay product 
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is assumed to be a neutrino of zero rest mass. The energy distribution 
given by Michel (1950) can be used to obtain the momentum distribution, 


Pree ig at" (WJ (pF)? 
PD) P= TR one WV oH) EIRP 


x {BM (W—/(p?+ 2?) +p?) [2M (W—+/(p?-+p2))+322]+4 32} 
NG lel Vinge Tome eee aaa §) 


W=(M*—+-p?—)°)/2M is the maximum total energy received by the 
#-meson in the centre of mass system of the disintegration. 


Fig. 1 


No. of Events 


ce) Kee} 200 300 
P (Mevé) 


Momentum distribution of the charged secondary in the centre of mass system. 
(The solid line represents the case V+>yu+-+7°+v; the dashed line 


Vt +yu++v-+y.) 


This function is plotted as a solid line in fig. 1 for the reaction (1) given 
above. When w is taken as the p-meson rest mass and A that of a 
7°-meson, the upper momentum limit of 330 Mev/e corresponds to a 
V*-particle rest mass of 1380 m,. This upper momentum limit has been 
chosen to give a reasonable fit with the data. 

For the decay scheme (2) we set A equal to zero. Equation (4) then 
can be simplified to give 


DA 


F(p) dp= a [BMP +p?)—6M v/(p? +p") + 2p" |p" dp. (5) 
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This function is plotted as a- dotted line in fig. 1. As before tug upper 
limit has been taken to be 330 Mev/c. This corresponds to a V -particle 
rest mass of 1330 m,. 

It is clear from fig. 1 that the momentum spectra for the two cases 
are very nearly the same except at the upper momentum limit. At the 
limit the dashed curve has a finite value and the solid curve goes to zero. 
Because the two curves have such similar shapes, one can be sure that 
when they are used to average the angular distribution, very similar 
results will be obtained. Hence only the distribution given in eqn. (5) 
will be used to calculate the desired transverse momentum distribution. 


(c) The Transverse Momentum Distribution 


The transverse momentum distribution is given by the expression 
Pmax pr 


= oe dp. <t ) a A 
Pon | arpa Pe) (6) 


Using eqn. (5) above for F(p) and integrating we obtain the result 
Ppp) dpr=K {V(x Ps") 3AP-+42)— 3 +8 Phas + Pa?) 


BLUE W+V(Pinax—Pr") 
— 24142 Sooo, een 
5} (pr +H") In W—V/(Pax— Pr’) Pr app, (7) 


where K is a constant, and W= 4/(p>ax+p). 

The expression P(p7) is plotted in fig. 2. The area under the curve 
is normalized to be the same as the area under the histogram of the 
combined experimental data given by Armenteros et al. (1952) and 
Barker et al. (1952). 


§ 3. CONCLUSIONS 


(a) The V+-particle data of Armenteros ef al. (1952) is consistent with 
a single three-particle decay scheme. 

(b) The transverse momentum distribution of the charged secondaries 
is insensitive to the assumed decay scheme of the V*-particle. 

In view of the very small amount of data available and the difficulty of 
rapidly obtaining more, a somewhat different line of approach seems 
necessary to distinguish between the modes of decay. Michel (1952) 


has pointed out that if the decay scheme (2) above actually occurs, then 
one would expect the reaction 


Ve => ety 2 Vethen! eee (8) 


to compete with it and have a comparable lifetime. However, no electron 
secondary has been reported thus far in the literature. This does not 
rule out (2) as the true mode of decay, but is an argument against it. 
On the other hand, the 7°-meson in scheme (1) has never been observed and 
an exceptionally clear case will be required to identify unambiguously an 
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electronic cascade resulting from the 7°-decay associated with a charged 
V-particle. Such an event or a clear case of an electronic secondary as 
indicated in (8) would be more useful for deciding which is the true scheme 
of disintegration than a study of the momentum spectra. 

The excellent agreement between eqn. (7) and the data of the Pic-du- 
Midi group, as shown in fig. 2, must be considered fortuitous. Because 
the histogram has been made up from only eighteen events, the possible 
statistical fluctuation of any point on the histogram is large. It should 
also be noted that if there is more than one type of charged V-particle, 
it would be possible to fit the observed transverse momentum distribution 


Fig. 2 


Na of Events 


° 100 200 300 
R (Mey) 


The transverse momentum distribution of the charged decay product of 
V+-particles which decay in flight. (The solid curve represents the decay 
scheme Vt>p+-+v-+-v and the histogram gives the V+-particle data 
of the Pic-du-Midi group.) 


by a suitable mixture of several two body or two and three body decay 
schemes. That such a mixture might be necessary is indicated by the 
recent announcement of O’Ceallaigh et al. (1952) that a new charged meson 
which decays into two particles has been discovered. It may well be that 
the charged V-particles are a mixture of these new mesons and the 


«-mesons found by O’Ceallaigh (1951). 
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Tr is generally well known that the present theoretical analysis of electron 
transport phenomena in metals rests upon an integral equation (the 
Boltzmann equation) for the electron distribution function. To solve this 
equation various assumptions and approximations, sometimes relatively 
crude, have to be made which we will not enumerate here in detail. As 
regards the thermal vibration of the lattice, however, the metal is assumed 
isotropic (when we may speak of longitudinal and transverse modes of 
oscillation) and in general the lattice wave velocity is treated as a constant. 
‘These two approximations indeed constitute the Debye model for the 
spectrum. 

Based on these assumptions previous treatments, e.g. Wilson (1937), 
Sondheimer (1950), have always predicted a minimum in the thermal 
conductivity of a pure monovalent metal at about 77=0-20. The experi- 
ments of Berman and MacDonald (1951) on very pure sodium showed no 
trace of such a minimum. Seeking the cause of this discrepancy we are 
examining carefully the assumptions used in the theory. 

As a first step, still taking the metal to be isotropic the transport 
equation has been solved without taking the lattice wave velocity to be 
constant (ef. Cornish and MacDonald 1951). For the thermal conductivity 
our calculations show that a minimum is still predicted. Hence the con- 
stant wave velocity assumption is not in itself the cause of this discrepancy 
between theory and experiment. 

We hope to publish in greater detail later. 
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The Mean Knergy of Three-Track Nuclear Disintegrations Produced by 
Cosmic-ray Particles 
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In 1949, Brown, Camerini, Fowler, Heitler, King and Powell measured 
the total energy of the particles emitted from disintegrations caused 
by cosmic rays as a function of star size for stars of from 4 to 18 tracks. 
They fitted their results with the empirical formula H=4N?+-37N, 
where FE is the total energy emitted in low- or medium-energy particles, 
which, if charged, make black or grey tracks (i.e. those with ionization 
greater than 1-5 times the minimum ionization of singly-charged particles), 
and N is the number of such charged particles. However, these 
investigators missed a considerable number of the small stars (Barford 
and Davis 1952). The stars missed are likely to be those of higher energy, 
which would considerably affect the above relation for small values of V. 
It therefore seemed worthwhile to repeat the measurements of the 
mean energy for small stars. 

Three-track stars in the G5 emulsion exposed under no absorber and 
-examined with high scanning efficiency by Barford and Davis were 
measured and their mean energy determined in the following way. 

The tracks of all particles of energy greater than 35 Mev which make 
an angle of less than 34° with the plane of the emulsion were grain 
counted, and the energies of the particles determined. Assuming that 
they were emitted isotropically the total number of such particles was 
calculated. As it was not possible to grain count tracks of particles 
between 25 and 35 Mev, their number was determined by fitting the 
energy spectrum of emitted protons determined by Camerini et al. (1950) 
to the observed distribution of particles with energies greater than 
35 Mev, and extrapolating back to the region between 35 and 25 Mev. 
The mean energy and number of evaporation particles less than 25 Mev 
was obtained from previous work using C2 emulsion (Harding, 
Lattimore and Perkins 1949, Hodgson 1952). The number of neutrons 
produced in these disintegrations was taken as 1-2 times the number 
of protons. The binding energies of the emitted particles were added 
to their kinetic energies and a correction made for deuteron and triton 
emission (Fowler 1950). 

As a result of this analysis the mean energy of the three-track stars 
was found to be 196-18 Mev. This may be compared with the value 
of 147 Mev obtained for three-track stars by extrapolating the formula 
of Brown et al. The difference between these two values may indicate 


that they tended to miss the more energetic disintegrations. 
$$ Ee eee 
* Communicated by Sir George P. Thomson, F.R..S. 
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A Note on the Measurement of Lifetime of Unstable Particles 
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Up to the present, only very approximate estimates of the lifetimes of 
the rare unstable particles, identified in recent years in cosmic ray 
phenomena, have been made. Several novel features are encountered 
when the problem of making a lifetime determination is examined. If 
the lifetime is longer than about 10~1° sec. there will be the possibility that 
charged unstable particles will come to rest in matter, and direct timing 
methods may be possible, provided that sufficient examples of decay can 
be observed and that the events under examination are not obscured by 
other phenomena. No comparable mechanism, however, is known which 
will bring uncharged particles to rest, and for these there seems little 
prospect of the application of direct timing methods. Up to now, it has 
not proved practicable to use conventional timing methods for any of the 
new particles, and the only data available which include a direct time 
element are cloud chamber photographs in which the act of decay is seen. 
The photographs also enable many of the different decay phenomena 
to be distinguished, and for most of the particles it is from these photo- 
graphs that measurements of lifetime must now be made. 

One kind of determination makes use of either measurements or estimates 
of the total flux of the particles under consideration, and the useful 
preliminary values of lifetime have in fact been reached in this way. 
We think that this approach is capable of only limited development, 
because it is so difficult to identify those particles which do not decay, 
and that, more definite information must be based explicitly on the data 
contained in actual records of identifiable decay processes. 


* Communicated by the Authors. 
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An examination of typical data—for example that of Armenteros ef al. 
(1952) and of Astbury et al. (1952)—shows that single groups of workers 
collect suitable data slowly, and so it is important that it shall be possible 
to combine data from many groups. In this note we suggest the way in 
which observations should be reported in order that they may be combined 
with others in a straightforward manner. 

Each photograph in which the decay of a particle is recorded is of the 
following type: at some point the unstable particle entered the useful 
volume of the cloud chamber (i.e. the part of the chamber which is suffi- 
ciently illuminated and which is within the depth of focus of the recording 
camera), and at a second point, if it had not decayed, it would have left this 
volume. The distance which it would have traversed between these two 
points, 2), may be described as the ‘ potential path length’. In fact, in 
each photograph decay took place after a ‘path length’, x, which is 
shorter than %. The data required from each event consists of the lengths 
X, %, together with additional information which will allow these lengths 
to be transformed to time intervals ¢, 4. Uncertainties of measurement 
are particularly serious for the transformation factor and for this reason we 
propose that the acceptable form of reporting observations should give 
separately the path lengths x, 2) and the observed features upon which 
the transformation factor may be based. 

It must be noted that the distances x, 2 cannot be measured to the 
extreme limits of the useful part of the chamber, for some short distance 
is required within which the act of decay may be recognized. These 
distances must therefore be measured inside fiducial surfaces located at 
some predetermined distance within the useful limits of the chamber, and 
reported measurements should give clear indication what reduction, if any, 
of this kind has already been made. For charged particles (V*, 7), the 
direction of motion of the primary particle is well defined, but for un- 
charged particles (V°, ,) this is not so, and it then becomes necessary to 
state the information upon which the trajectory of the particle has been 
determined. This may be based, for example, upon the identification of 
the source from which the particle has come, or upon the momenta and 
geometry of the secondary particles, together with the assumption of a 
particular decay scheme. 

The transformation from the measured distances 2, X) to the proper time 
intervals t, fj may occasionally be made directly when the ionization of a 
(charged) primary allows an estimate of the velocity of the particle to be 
made (see Astbury e¢ al. 1952), but more often it will be derived from the 
measured momentum of the primary and an assumed value of its mass, 
or estimated from measurements made entirely upon the secondary 
particles. The reliability of the transformation factor will certainly 
vary widely from one example to another, and it is particularly important 
rele aa Mees basis of each value should be reported adequately. 
are a ee et ene of the traversal times f, ty the most probable life- 

ced following the procedure discussed by Peierls (1935) 
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and Hole (1947). A statistical treatment has been given by Bartlett 
(1936). It is readily shown that if the average value of a/2, is close to one 
half, then the lifetime, 7, is probably longer than the average value of ty, 
and a very large amount of data will be required before a reasonably 
accurate value of 7 can be deduced. It will be shown in later publications 
that this is the situation which obtains for V*-particles, for some at least 
of the apparatus now in use. For these particles, larger values of ty 
are accordingly desirable, and may be obtained either by using a larger 
chamber or by modifications of selection which increase the number of 
slow particles recorded. 
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Hyperfine Structure in Neodymium Ethyl Sulphate 


By B. Buranezy, H. E. D. Scovim and R. 8. TRENAM 
Clarendon Laboratory, Oxford * 


[Received July 9, 1952] 


HYPERFINE structure in the paramagnetic resonance spectrum of neo- 
dymium ethyl sulphate was first observed by Bleaney and Scovil (1950), 
who analyzed the spectrum in terms of the usual Hamiltonian (with 
S=4, [=7/2 for the two odd isotopes 143 and 145). 


7-09, 11,8, pg (1 8,-- H,S,)+AS,1,+ B(S,1, 8,1) 
ar ee oe SAC 1)}. 


Further measurements of the strong field spectrum at 3 cm wavelength, 
using proton resonance to measure the magnetic field, have yielded the 
following values for the parameters (A, B in units of 10-4 cm~*) : 


gu=3:535-+40-001 — g, =2-073-40-002 


Isotope 143 Isotope 145 Ratio 
| A] =380-3+0-1 A =236-4-+0-1 1-6083 +-0-0012 
| B| =198-9-+0-5 B=123-7 0-5 


The values of B are less accurate, partly because not all the hyperfine 
structure is fully resolved, and partly because the theory has not been 
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carried far enough. The energy levels, with the magnetic field perpendicular 
to the axis, have been evaluated by fourth order perturbation theory. 
Asmall discrepancy between theory and experiment remains in the spacing 
of the hyperfine structure, which might be due either to higher order 
effects or to a nuclear quadrupole term with | P| of the order of 2-10~4 em™*. 
The latter possibility has been eliminated by making observations at 
3cm wavelength of weak transitions corresponding to changes dm= = 1 
in the nuclear magnetic quantum number, which should be split into 
doublets by the nuclear quadrupole interaction (Bleaney 1951). No 
such splitting could be detected, giving an upper limit for P for either 
isotope of 1 10-4 cm}. 

To obtain further information about P the hyperfine transitions have 
been plotted into zero magnetic field using wavelengths of from 7 to 
20cm. At zero field, eight strong transitions are allowed for each 
isotope, three of which involve only the two parameters (A—2P)? and 
B?. By taking these transitions in pairs, three values of each of these 
parameters are obtained which are just consistent within the experi- 
mental error. They are (in units of 10-4 cm~?) 


(A 2P | | B| 
143 145 143 145 
383-0-+41-0 233-7+1-0 197-6-+1-0 125-5-+1-0 
381-3-+41-3 235-2+1°5 198-4-+0-4 124-6 0-5 
382-5 40-6 234-0-+0-6 198-5 40-4 124-7-+0-5 


By combining the weighted average of these values with measurements 
of transitions involving only A? and B?, the following values for all the 
constants were obtained. 


143 145 Ratio (143 : 145) 
|A| =381-3+0-3 | A| =235-5+0-2 1-619 0-003 
| B| =198-5+0°3 | B| =124-9+0-5 1-590 +0-007 
|P|= 0-7-0:3 [P| = 0-8+0:3 


It will be seen that the zero field measurements are inconsistent in 
that the ratios of A and B differ from each other and from the strong 
field value. The sign of P, relative to A, is opposite for the two isotopes, 
but little significance can be attached to this, or to the actual magnitude 
of P, until the reason for the discrepancies has been established. At 
present no satisfactory explanation has been found, but it is planned to 
repeat the zero field measurements with a crystal grown from heavy 
water. The reduced line width obtained thereby should increase the 
accuracy sufficiently to provide a good test of the spin Hamiltonian used 
to represent the results. 

These small inconsistencies do not affect the size of the nuclear magnetic 
moments, in which the main uncertainty lies in the value of fr. 
Taking 1/r=34-10%4 cem-3, Elliott and Stevens (1950) give p44,;=1-4, 
/145—0°85 nuclear magnetons. It seems probable that the actual values 
of P are not greater than 1-104 em-}, corresponding to an upper limit 
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for the nuclear quadrupole moments of 1-2 10-24cm2?. The signs of the 
nuclear magnetic moments are not determined by the measurements; 
but from the optical spectrum of the separated isotopes, Murakawa and 
Ross (1951) give 443 —1-0-+0-2, 441 4,= —0-62-++0-09 nuclear magnetons. 
The discrepancy between these values and ours may be due partly to 
the fact that, in the optical case, all contributions to the hyperfine 
structure except those from the 6s electrons are ignored, and partly to 


uncertainty in the value of 1/r°. 
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On the Production of «-mesons 


By U. Haser-Scuaim and G. YEKUTIELI 
The Weizmann Institute of Science, Rehovoth, Israel* 


[Received June 22, 1952] 


RECENT experiments of the Bristol group (Daniel ef al. 1951) have 
indicated that heavy mesons («-mesons or kappons, for short) are found 
in appreciable numbers. In some stars produced by very energetic 
primaries (>50 Bev) a ratio between the number of kappons and 
pions n,/n,—=0-6 was claimed. Such a large ratio suggests that the 
kappons like the pions are strongly coupled to the nucleons. 

One may try to account for the production of the kappons by either of 
the following assumptions : 

(a) The kappon is a Yukawa type particle which interacts with the 
nucleon through a ‘kapponic charge’ quite independently of the 
existence of the ‘ pionic charge’. This would lead us to believe for 
example, that the force between two nucleons at rest is derivable from 
two Yukawa potentials with ranges #/m,c and h/m,c. 

(b) The production of both kinds of mesons is the result of one field. 
In this case neither of the mesons is a Yukawa particle in the full sense 
of the word. 

We shall not investigate these assumptions here but rather limit 
ourselves to their manifestation in the case of very strong coupling when 
Fermi’s statistical theory is applicable (Fermi 1950, 1951). 

In this extreme case assumption (a) will lead us to the conclusion 
that the kappons are produced in a volume of the order of 


Ag [f heX3 
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whereas the volume available for pion production will be 


Hence even at the highest energies (7/7) <[(m,/m,)*] except for possible 
spin factors. On the other hand, assumption (b) will lead to 2Q,=2,. 
In this case the high energy limit of n,/n, will be determined by spin 
factors only. 

The number of kappons ,, and pions 7,, produced in nucleon-nucleon 
collisions is calculated in the following way. Both kappons and pions 
are taken to be bosons with spin zero which can be either charged or 
neutral. Only the pions are treated as extreme relativistic particles in 
c.m, system whereas for the kappons the exact relation between energy 
and momentum is used and the integrals in momentum space are 
computed numerically. The results of the calculation are summarized 
in table 1, where W’ is the energy in the laboratory system. The following 
values for the masses are used: m,=0-69 and m,=0-15 in nucleon 
mass units. 


Table 1. n, and», as Functions of Energy of Incident Nucleon 


W’ in Assumption (a) Assumption (b) 
Mc? 

Nx Nx Ny | Nn Nx Nx Ny/| Nae 
20 0-011 2-60 0-0042 0-63 1-78 0-35 
50 0-019 3°63 0-0052 Aly 2-34 0-50 
100 0-027 4-55 0-0060 1-64 2-85 0-57 
200 0-038 5:53 0-0069 2°25 3°43 0:66 
300 0-046 6-20 0-0074 2-64 3°85 0-69 
500 0-057 7-18 0-0080 3-17 4-38 0:72 
700 0-065 7:80 0-0084 3°56 4°75 0-75 
1000 0-075 8-50 0-0088 4-05 5-10 0-79 


Table 1 clearly favours assumption (6). 

Additional evidence concerning the two assumptions can be drawn 
by comparing the interaction mean free paths of kappons, /,, and pions, J, _ 
in photographic emulsions (compare Daniel et al., p. 22). Under 
assumption (a) one obtains a ratio (i,/l,)=5-4 whereas (I,/l,)=1 under 
assumption (0). 


REFERENCES 


DANIEL, R. R., Davies, J. H., Mutvery, J. H., and Perxrs, D. H., 19513 
Bristol Conference of V Particles and Heavy Mesons. 
Frrmt, E., 1950, Prog. Theo. Phys., 5, 570 ; 1951, Phys. Rev., 81, 683. 


Correspondence 999 


Paramagnetic Resonance in Praseodymium Ethylsulphate 


By B. Bueanzy and H. E. D. Scovm 
Clarendon Laboratory, Oxford* 


[Received July 21, 1952] 


THE paramagnetic resonance spectrum of praseodymium ethylsulphate 
in a crystal diluted with lanthanum has been observed at 20° x, using 
radiation of wave numbers 0-332, 0-687 and 0-789cm-1. A single 
electronic transition is observed, with a hyperfine structure of six equally- 
spaced components. This confirms the known nuclear spin (5/2) of the 
only stable isotope, 141. 

The ground state of the free Pr+++ ion is 4f?, 3H,, which, under the 
action of a crystalline electric field of C;, symmetry, splits into a singlet 
and four doublets. The theory of Elliott and Stevens (unpublished) 
shows that a doublet, whose states may be written as 

cos 6 |J,=+4)+sin 6| J,=F2) 
(the upper sign refers to one state, the lower to the other), should lie 
at the bottom. This is not a Kramers’ doublet, and, since the states 
contain no values of J, differing by unity, the normal type of resonance 
transition is not allowed. If there is any dissymmetry in the crystalline 
field which admixes the two states, however, transitions are allowed with 
the r.f. magnetic field parallel to the crystalline axis. We have tentatively 
tried to allow for such an effect by adding a term in S,, to the usual spin 
Hamiltonian. It would then be written 

CE Jnplts A ee, ee a ee oo CL) 

with an effective spin S=}. The absence of matrix elements in the mag- 
netic field between the two states of the doublet is indicated by the 
absence of the usual terms g,B(H,S,+H,S,) and B(S,J,+S,1,). This 
Hamiltonian gives energy levels 

W=1t}(4?4 {9,8o,+Am¥), 2 2. 2 2. . (2) 
where m is the nuclear magnetic quantum number, and transitions are 
allowed only between states of the same m. (We have omitted the nuclear 
electric quadrupole term, as this is likely to be negligible compared with 
the other terms.) 

Observations at three wave-lengths show six lines, equally spaced in 
magnetic field, as would be expected from the form of eqn. (2), when m 
is given the values 5/2, 3/2,....—5/2, W is kept constant, and H is 
varied. The field corresponding to the centre of the spectrum is not 
quite proportional to the frequency of observation. This indicates that 4 
is very small (approximately 0-04cm™1) and gives g,=—1-69--0-01, 
A=0-083-+0-001 cm-1. When the axis of the crystal is rotated away 
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from the magnetic field, the spectrum moves rapidly into higher fields 
and out of the range of our magnet, indicating that g, is very smal] 
(<0-3). Susceptibility measurements of van den Handel (1941) indicate 
also that g,=0, but give a slightly smaller value (1-56) for g,. 

This discrepancy may be due to one of two causes : 

(a) The resonance lines are several hundred gauss broad and asym- 
metrical in shape, the intensity rising slowly on the low field side and 
falling steeply on the high field side. We have assumed that the point of 
maximum intensity gives the position of the unbroadened line. The 
asymmetrical shape would, qualitatively, be expected from the form of 
eqn. (2), but only if 4 were comparable with the frequency of observation. 

(b) Our assumed Hamiltonian (1) may be wrong. The admixture of 
the two states and consequent splitting of the electronic doublet is 
presumably associated with the Jahn-Teller effect, which lifts the 
degeneracy by allowing a distortion of the crystal field which lowers the 
total energy of the system. It is doubtful whether the addition of a simple 
term such as 4S, fully represents this complex process, which is little 
understood (cf. Bleaney and Bowers 1952). 
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The Dispersion Formulae and the Polarization Fields 


By Sir K. 8. Krisunan, F.R.S., and 8. K. Roy 
The National Physical Laboratory of India, New Delhi* 


[Received June 18, 1952] 


In the usual derivation of the Drude formula for the dielectric constant 
of a dense medium as a function of the frequency of the incident waves, 
since one does not invoke the presence of a polarization field, and since 
further it is known that the Lorentz dispersion formula reduces to the 
Drude formula in the special case when the polarization field is made to 
vanish, any verification of the Drude dispersion formula is. sometimes 
taken to imply the absence of a polarization field.» This conclusion is 
not justified, and it is the main purpose of this note to emphasize that 
the Drude formula for the dielectric constant as a function of the frequency 
of the incident waves (as distinguished from the formula for the dielectric 
constant as a function of the density) is perfectly consistent with the 
occurrence of a polarization field having the Lorentz value, or any other 
value. This will be the case even when the polarization field factors 
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pj; defining the polarization field acting on any oscillator of type 7 due 
to all the oscillators of type j in the medium differ from the Lorentz 
value 47/3, and differ from one another. In other words, irrespective 
of the actual polarization field that occurs in the medium, i.e. irrespective 
of the actual values of p,, that define this field, the observational data 
for the dispersion of the dielectric constant can always be represented 
by the Drude formula 
A. 
ae ae LY RCT SE Nein a em seen 

the characteristic frequencies w; that appear in the formula being the 
resonance frequencies of the medium. 

The same observational data can also be fitted alternatively into a 
formula of the type 


Cee = AGN eT Ve eS Te ae Go, os, ods y(t) 
where ; 
B; 
Ki= Diet ePuxy), + + + +e = (8) 
1 j 
in which j can take all values i,j, k,.... When all the p,,’s have the same 
value, equal to p say, Xpyx;=px, and eqn. (2) reduces to the form 
3 
qr C; 
eee OFA? (4) 
where 
4a 
Ci a Re A oe ee G3) 
, (5) 
and 
C= pB,. (6) 


When p has the Lorentz value 47/3, it will be seen that « reduces to 2 
and eqn. (4) to the simple Lorentz formula. For convenience we shall 
describe the type of dispersion formula defined by (2) and (3) as a 
generalized Lorentz formula, as distinguished from the simple Lorentz 
formula corresponding to «=2 in eqn. (4). 

Thus any observational data on the dispersion of the dielectric constant 
of a dense medium can be fitted into a Drude formula, or alternatively 
into a generalized Lorentz formula involving any given set of polarization 
factors p;, This fitting into the alternative formulae can be done 
irrespective of what the actual polarization field in the medium may be, 
and even if there is no polarization field at all. 

The equivalence of the Drude and the simple Lorentz formulae for the 
dispersion of the dielectric constant of a dense medium was first 
demonstrated by Livens (1912) for the case when the medium has a 
single resonance frequency, and by Herzfeld and Wolf (1925) for the case 
when the medium may have more than one resonance frequency. The 
equivalence extends, as we have just seen, also to the generalized Lorentz 
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formula corresponding to any specified set of polarization factors Pij 
defining the interactions between the dipole moments of the different 
oscillators. 1 

This equivalence arises from the following circumstance. The effect 
of any polarization field that may be present in a dense medium on its 
dielectric constant may be taken into account in two alternative ways, 
which are equivalent. The effective field that polarizes the medium 
may be taken to include, in addition to what is usually defined as ‘ the 
field in the medium’, the polarization field also, in which case the 
dispersion formula will be of the Lorentz type and the characteristic 
frequencies that appear in the formula will be the frequencies {2; of the 
individual oscillators, regarded as unaffected by the presence of the 
polarization field, ie. unaffected by their mutual interactions. Alter- 
natively one may also regard the polarization field as effective in changing 
the frequencies 2; to the corresponding resonance frequencies w,; of the 
medium, and its effect on the dielectric constant as exercised indirectly 
through these frequencies. The result is the Drude dispersion formula. 
(The important part played by the polarization field in determining the 
infra-red resonance frequency of the alkali halide crystal, in which the 
polarization field involved is readily calculated, is discussed by us in a 
recent paper (Krishnan and Roy 1951).) 

Conversely, starting with the Drude formula, and the known resonance 
frequencies w, that appear in the formula, and taking them to correspond 
to zero polarization field, one may postulate any desired polarization 
field and reduce the Drude formula to one of the Lorentz type. The 
frequencies 2; that appear in the latter formula will be different from 
the corresponding resonance frequencies w; by amounts that will be 
determined by the polarization field postulated. 

Hence dispersion data as such can not give us any information regarding 
the nature of the polarization field that exists in the medium, though 
they enable us to obtain the resonance frequencies, directly if the data 
are expressed in the Drude form, since the characteristic frequencies 
w,; that appear in this formula are just these frequencies, or after a 
simple calculation from the corresponding Q; if the data are expressed 
in the Lorentz form. 

The position, however, is very different with expressions for the 
dielectric constant as a function of the density, where the equivalence 
of the Drude and the Lorentz types does not hold. 
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The X-ray Spectrum from a 70 MeV Synchrotron 


By I. B. McDrarmip 
Synchrotron Laboratory, Queen’s University, Kingston, Ontario, Canada* 


[Received July 11, 1952] 


THE spectrum of the x-rays produced in the forward direction by 
70 Mev electrons striking a conventionally thick target is being measured. 
The scattering method has been used to determine the energies of 
637 electron pairs produced by the x-rays in G5 emulsions and from 
this the spectrum has been deduced. The number of photons found 
below 25mMev is considerably less than that predicted by the Bethe- 
Heitler—Schiff (1951) theory corrected for target thickness. 

The target in the synchrotron here is a tungsten rod about 4in. in 
diameter. Presumably the x-rays are produced by what one might call 
glancing blows of the electrons on the target, and hence the effective 
thickness is not known accurately. However, measurements by Holloway 
and Janzen of the angular width of the beam indicate a target thickness 
of the order of 0-05 to 0-1 cm (Schiff 1946). The effect on the spectrum 
of a thick target has been calculated and a result, similar to that given 
by Eyges (1951), obtained. The main effect of the calculation is a 
relative lowering of the high energy end. The best agreement with 
theory, above 25 Mev, is obtained if we assume a 0-1 cm thick target. 
The smooth curve of fig. 1 is Schiff’s integrated photon distribution 
corrected for this target thickness. 

The G5 emulsions were placed in the centre of the beam 3 m from 
the target and given approximately 10 milli-r of radiation. After 
processing the plates were searched systematically for electron pairs and 
from the measured mean angle of scattering the energies calculated. 
The method used to measure the scattering and the evaluation of the 
constants in Williams’ (1940) scattering theory for G5 emulsions have 
been described previously (McDiarmid 1951). 

The pairs used for the spectrum determination were selected in the 
following way. The scattering was measured on all pairs in which both 
electrons travelled at least 300, in the emulsion. Those pairs were 
discarded for which the probable error in the energy determination of 
either electron (obtained from the measured ‘ noise level’ and the length 
of track in the emulsion) would have been greater than 15%. This 
selection process discriminates slightly against the lower energy pairs 
and this was corrected for as follows. For a given energy and the 
measured ‘ noise level’ the length of track required to give a probable 
error in energy less than 15% was calculated. From the distribution 
function of Rossi and Greisen (1941) and the known emulsion thickness 
Oe aE Sl oh Eo le ee a Ne al Bc oe ee 
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the fraction of electrons that escape before going this distance was 
calculated and from this the corrections to be applied to the number of 
pairs in the various intervals obtained. The total corrected number 
of pairs was approximately equal to the total number found in the 
emulsion indicating that the corrections applied were of the right order. 


Fig. | 


THEORETICAL PHOTON DISTRIBUTION 


NUMBER OF PHOTONS (ARBITRARY UNITS) 
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PHOTON ENERGY (MEV) 


Number of photons as a function of energy. The theoretical curve is Schiff’s 
integrated distribution corrected for a 0-1 cm thick target. 


The corrected number of pairs in the different energy intervals was 
divided by the corresponding cross section for pair production, giving 
numbers proportional to the number of photons in the beam at these 
energies. The measured photon distribution obtained in this way is 
compared with the theoretical distribution in fig. 1. Evidently, there is _ 
a large discrepancy between the measurements and theory below 25 Mev. 

The photon spectrum we have found is very nearly independent of 
energy in the region 15 to 25 Mev and this is in agreement with ‘ relative 
yields ’, as defined by Perlman and Friedlander (1948), of (y—n) reactions 
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measured in this laboratory by Edwards and MacMillan (1952). The 
definition of yield leads to the approximate expression 

Y, NA, 

PeemN emesis MY) 
for the ratio of two yields, where Y is the yield measured here, N is the 
number of photons per Mev around the average energy producing the 
reaction, and J is the integrated cross section. From the yields measured 
in PC, Cu, Cu, and Ag the ratios V ,/N, for several pairs of energies 
have been calculated from eqn. (1). The average energies producing 
the reactions and the integrated cross sections were taken from the work 
of Katz and Cameron (1951), and Diven and Almy (1950). In table 1 


Table 1. Ratios of the Number of Photons Around Different 
Pairs of Energies. 


No. at 18 Mev} No. at 17 Mev| No. at 15 Mev| No. at 15 Mev 
No. at 24 Mev] No. at 24 mev] No. at 24 mev| No. at 18 Mev 


From yield 
measurements 0:97 1:0 1-0 1:0 


From pair 
measurements ~l ~l ~l ~l 


From Schiff’s 
theory 


these ratios are compared with the corresponding ratios given by the 
pair measurements and by Schiff’s theory. It can be seen that the 
yield and pair measurements predict about the same ratios while the 
theory predicts much higher ratios. 

A complete report of this work is being submitted to the Canadian 
Journal of Physics. 


I wish to thank Professor J.-A. Gray for suggesting this problem and 
for many helpful discussions. I am also indebted to the National 
Research Council of Canada for financial assistance. 
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Interpretation of the de Haas—van Alphen Effect 


By L. ONsSAGER* 
The Royal Society Mond Laboratory, University of Cambridget 


[Received July 14, 1952] 


Ar low temperatures many metals exhibit an oscillatory variation of 
their magnetic susceptibilities with the field intensity. This effect was 
first discovered by de Haas and van Alphen in a study of bismuth. 
Recent investigations have shown that it is quite common ; the remaining 
negative findings might well be due to the limitations of experimental 
technique. 

The effect is generally attributed to the helical motion of the free 
electrons, which ought to be quantized as far as the components of 
momentum perpendicular to the magnetic field are concerned; the 
component parallel to the field remains a constant of motion capable of 
all values. The theory applicable to free electrons was elaborated by 
Blackman (1938), Shcenberg and Landau (1940), who recognized the 
effect of the periodic field in the crystal by the introduction of modified, 
generally anisotropic, effective electron masses. According to these 
theories oscillations of long period, which are studied with relative ease, 
can occur when two of the principal masses are small. The zone theory 
permits this in rather special locations in the reciprocal lattice, most 
commonly perhaps where Bragg reflections of equal strength intersect 
(Robinson 1950), but hardly over extended domains in momentum 
space ; wherever a small effective mass is found, its constancy is in doubt. 

Under the circumstances, a more general approach to the problem is 
of interest. We retain the wave-packet picture, which is more or less 
implicit in previous theories. Let the energy in the absence of external 
electric and magnetic fields be a function W(p) of the momentum p ; 
then the group velocity v is given by the relation SW=v. dp, and an 
external force X will cause an acceleration given by p=X. 

Admittedly, the effect of a magnetic field cannot be described quite 
rigorously in these terms. We replace p by the ‘kinetic momentum’ II, 
which contains the vector potential A: 


T=p-+(c/c)A 


but the statement W—W/(II) is somewhat ambiguous because the 
components of II do not commute. It would appear that a comprehensive 
analysis of the difficulties which arise thus has not yet been achieved. 
In kinetic theories (Shockley 1950) it is customary to ignore such 
complications, and it does seem reasonable to hope that neither previous 
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theories of diamagnetism nor the present generalization will be invalidated 
by the errors involved in the ‘ wave-packet ’ approximation, or at least 
that the error in the computed susceptibility will not vary rapidly with 
the field intensity. 

According to our picture, then, a magnetic field causes the vector II 
to move along the line of intersection between a plane (II . H)=constant 
and a surface W(II)=constant in momentum space (Shockley 1950). 
The helical path of the electron, projected upon a plane perpendicular 
to the field, becomes an exactly similar curve in ordinary space, turned 
through a right angle and of dimensions changed in the ratio c/eH. The 
periodic component of the motion, which involves the components of 
p and r perpendicular to H, is quantized. We apply the Bohr-Sommerfeld 


tule in the form 
bp .dr=(n+O)h. 


Then the cross section of the helical orbit is determined by the condition 
that the enclosed magnetic flux equal (n+ @) ‘ force lines’, whereby we 
recognize a fundamental unit 


1 force lne=hc/e=4rg=4-135 x 10-7 maxwell 


equal to the flux from one of Dirac’s hypothetical magnetic poles. The 
energy levels are degenerate because the location of the helical axis 
is arbitrary ; the multiplicity of each level equals the number HQ/4zg 
of force lines traversing the cross section 2 of the specimen. The area 
enclosed by an ‘ orbit ’ in J/-space is proportional to the field intensity, 


(n+0)(he/c)H ; (momentum units) 
or (n+ @0)(e/hc)H ; (wave-number units). 


The known explanation of the de Haas—van Alphen effect depends on 
this, that with increasing field intensity, the electrons are redistributed 
among a decreasing number of levels. The capacity (degeneracy) of 
each level increases in proportion with the field intensity : the spacing of 
the levels becomes correspondingly more open. The top levels become 
depleted in succession, and the magnetic susceptibility will vary in 
phase with the state of depletion of the highest level which retains an 
appreciable number of electrons. In the cases which have been accessible 
to study by known techniques, we have reason to believe that the top 
energy level £ of the electrons was kept sensibly constant by reservoirs 
of electrons which would not produce appreciable oscillations at the 
prevailing field intensities. Under those conditions, when the intensity 
of the field is increased, the nth quantum state will be (more than half) 
depleted as soon as its energy exceeds ¢. The field necessary to 
accomplish this depends simply on the area in momentum space enclosed 
by the curve W(p)= on some plane (p.H)=constant. If the area in 
question does not depend on the component of momentum parallel to H, 
that is, if the surface W(p)=¢ is a cylinder, then the depletion will take 
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place simultaneously along the entire length of the cylinder and cause 
rather large oscillations. This improbable contingency does not seem 
to be realized in any case which has been studied. If the Fermi surface 
W(p)=¢ is convex, then depletion will start slowly at the ends and 
proceed at an accelerated pace towards the middle. The rapid last stage . 
of depletion will correspond to the maximum area of intersection, and 
in line with the known theoretical results for ellipsoidal surfaces, we 
shall observe oscillations whose period depend on this maximum area. 
Alternatively, the Fermi surface may resemble a hyperboloid of one 
sheet, in which case we shall expect oscillations which correspond to the 
minimum area of intersection. It is even conceivable that the Fermi 
surface may take the shape of an endless tube of variable cross section 
spanning the entire length of a cell in wave-number space. In that case 
we have to expect superposed oscillations of two different periods, which 
will give us measures of the smallest and the greatest cross section, 
respectively. According to preliminary results obtained in this 
laboratory, such a situation may prevail in graphite. 

.To summarize: the oscillatory diamagnetism reflects specific 
geometrical features of the Fermi surface; the nth maximum (counted 
from 1/H=0) occurs for a field intensity H, given by the relation 


n+@=(he/e)A/H,,. 


where A is any extremal area of intersection between the Fermi surface 
and the family of planes (H . p)=constant. 

The phase 9 contains a known correction which depends on the signs 
of the principal curvatures of the Fermi surface. Otherwise, the present 
considerations add nothing much to the existing theory as regards the 
phases and the amplitudes of the oscillations. Indeed, the generalized 
theory hardly foresees results which could not be interpreted in terms of 
ellipsoidal Fermi surfaces except this, that one negative principal mass 
is permissible and indicates a hyperboloid. 

In order to explain the observed amplitudes and phases in full detail, 
the theory has to invoke effects which are not fully understood. This 
situation remains unchanged. 
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C. Notices of New Books and Periodicals received 


Lhe Principles of the Control and Stability of Aircraft. By Prof. W. J. Duncan, 


D.Sc., F.R.S. [Pp. XVI+384.] (Camb. Aeronautical Series ; C. U.P. 1952.) 
Price 40s. net. 


THE pace of aeronautical progress is such that there is little time for the 
consolidation of knowledge, and still less for the writing of text-books. The 
appearance of a first volume in a Cambridge Aeronautical Series is therefore 
very welcome. The volume under review is an excellent book, and fills a 
serious gap in the literature: no major work exclusively devoted to stability 
and control of aircraft has appeared for forty years. It contains chapters on 
the mechanics of flight, on the equations of motion and their solution, on 
longitudinal-symmetric and lateral-antisymmetric motion, on flap controls and 
on aerodynamic derivatives. A chapter is devoted to each of the important 
subjects of static stability and manoeuvrability, distortion, and compressibility 
‘effects. Two chapters, on spinning and on landing flaps, by Professor A. D. 
Young, are also included. The reviewer found the book extremely informative 
and very readable, if a little academic in parts; it is excellently set out and 
printed, and is remarkably free from errors. ceive 


Electronic and Ionic Impact Phenomena. By H.S. W. Masszry and E. H. 8S. 
Buruop. [Pp. 669.] (Oxford University Press.) Price 70s. net. 


Owine to the severe limitations of space in this Journal, an adequate review 
of a closely written volume of 669 pages is clearly not possible. Nor can the 
reviewer pretend to have read it all at this stage. But it is a work that can 
be most strongly recommended, and one in which the authors have undoubtedly 
succeeded in their aim to give a balanced view of the subject from both the 
theoretical and experimental standpoints. It is indeed full of informa- 
tion ; yet in those portions which the reviewer has tested, because of a special 
interest in them, the data are critically selected and presented. Collisions of 
electrons and ions with atoms, molecules, surfaces and ions of opposite sign, as 
well as their passage through gases, all have chapters devoted to them, and the 
volume will serve as a valuable source of reference in many fields of work. 


oN fi ad 


The Auger Effect. By E. H.S. Burnor. [Pp. 188.] (Cambridge Monographs 
on Physics.) Price 27s. 6d. 


In his introductory chapter the author describes and _ illustrates the 
cloud-chamber events which were interpreted by Auger as radiationless 
re-arrangements of the electrons of an excited atom. There too are discussed 
certain other radiationless transitions, the most important of which are internal 
conversion of y-rays and internal production of electron pairs. Under ‘ Theory 
of the Auger Effect ’ are listed the formulae necessary for the calculation of 
the transition rates for the different phenomena. The equivalence of the two 
treatments of the Auger Effect—as an internal absorption of radiation and 
as a direct interaction between electrons—would be more easily seen if the same 
normalization of the wave-functions had been used in the two cases. An 
exhaustive summary of the experimental methods used in studying the effect 
is given. The various roles which the radiative transition is required to play 
in the interpretation of x-ray spectra are examined in detail. Useful collections 
of curves and tables of calculated internal conversion coefficients are included 
in the chapters on internal conversion. The influence of the Auger effect on 
meson capture and its relation to atomic and molecular spectra are treated 
more briefly. G. E. L-W. 
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The Alchemists, Founders of Modern Chemistry. By F. SHeRwoop TaYtor. 
[Pp. 246.] (William Heinemann.) Price 12s. 6d. 


Tux results of Dr. Taylor’s 20 years of study of the history of alchemy are 
distilled into this fascinating book. Although, as the author says, alchemy 
has been described as the pursuit of error, its study is useful in giving an insight 
into the history of the evolution of advanced techniques, and into the characters 
of serious experimenters whose most lasting contribution to science was the 
demonstration that their modes of thought could lead to useful results only by 
accident. 

Dr. Taylor describes this book as a sort of interim report on his findings—for 
the study of the history of alchemy is involved, confusing and controversial. 
Not only is it difficult for us nowadays to appreciate the underlying theoretical 
ideas which made their efforts to transmute elements seem plausible, but the 
alchemists deliberately concealed their most important instructions and results 
in a tangle of obscure symbolical language. Dr. Taylor makes the point that 
the transmutations that the alchemists dreamed of are possible nowadays in 
the atomic pile, and that this has indeed had precisely the effect that the 
alchemists feared—the placing of gigantic power in the hands of those who may 
misuse it. 

The book is illustrated with many quotations from original sources, often . 
accompanied by the translation into modern terminology. There are also 
valuable chapters on the transition from alchemy to chemistry—a change which 
was accomplished less easily than is often represented—and on the relation 
between alchemy and science. E. M. 


Schwerkraft und Weltall. By P. Joxpan. In Die Wissenschaft, Vol. 107, 1952. 
(Braunschweig : Vieweg und Sohn.) [Pp. 207.] 


THE first part of this book gives an introduction to Riemannian geometry and 
general theory of relativity. A welcome feature is the emphasis on physical 
interpretation. The second part covers the field indicated by the sub-title 
‘Grundlagen der theoretischen Kosmologie’. It is based mainly on the 
author’s own ideas about treating the gravitational ‘ constant’ as a field 
function in the 5-dimensional formulation of the Maxwell—Einstein theory. 

ALS. 


Mathematics, Queen and Servant of Science. By E.T. Brix. (G. Bell & Sons 
Ltd., London 1952.) [Pp. 437.] Price 21s. 


Tis book originated from two earlier popular books of the author’s ‘ The 
Queen of Science’, and ‘The Handmaiden of Science’. It presents very 
readable accounts of culminating points in past and present mathematics, with 
a special emphasis on mathematical logic, theory of numbers, algebra, the 
concept of invariance, and applications to mathematical physics. A.S 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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